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Abstract 



The non-Abelian analog of electrical conductivity at high temperature has 
previously been known only at leading logarithmic order: that is, neglecting 
effects suppressed only by an inverse logarithm of the gauge coupling. We 
calculate the first sub-leading correction. This has immediate application 
to improving, to next-to-leading log order, both effective theories of non- 
perturbative color dynamics, and calculations of the hot electroweak baryon 
number violation rate. 



I. INTRODUCTION 

We will provide a next-to-leading log order (NLLO) calculation of the non-Abelian (or 
"color") conductivity in hot, weakly-coupled non-Abelian plasmas. The motivation for this 
calculation, and an overview of the strategy, are presented in Ref. [1]. Here, we will simply 
get down to business. "Hot" plasma means hot enough (1) to be ultra-relativistic, (2) to 
ignore chemical potentials, (3) for non-Abelian gauge couplings to be small, and (4) to be 
in the high-temperature symmetric phase if there is a Higgs mechanism. 

As discussed in Ref. [1], there is a sequence of effective theories which describe color 
dynamics at large distance scales and long time scales. 

Theory 1: u,k < T 

(A + v-D)H/-v-E = 0, (1.1a) 
D V F^ V = f = m 2 (v^W) . (1.1b) 

These "hard-thermal-loop" equations amount to linearized, collisionless, non-Abelian, 
Boltzmann-Vlasov kinetic theory. The parameter m is the Debye screening mass, which 
is O(gT). This effective theory is valid, to leading order in the gauge coupling, for frequen- 
cies and momenta small compared to the temperature, u,k «T. 



Theory 2: u <C k <C m 

v-DW-v-E = -6CW + Z, (1.2a) 

(W)=0, (1.2b) 

DxB=j=m 2 (vW) , (1.2c) 

m = ^sc. (i.2d) 



This theory is a stochastic, collisional, linearized kinetic theory of hard excitations coupled 
to slowly varying non-Abelian gauge fields. Both the noise £, and the linearized collision 
operator 8C, arise from integrating out the effects of gauge field fluctuations below the scale 
of m. Here, and henceforth, ((• • •)) denotes an average over the (Gaussian) stochastic noise. 
Theory 2 is valid for spatial momenta small compared to the Debye screening mass and 
frequencies small compared to momenta, uo <C k <C m = O(gT). Eq. (1.2b) implements the 
effects of Gauss' Law in this range of cu and k and is explained in Ref. [2]. 

Theory 3: uj < k < 7 

(jE = DxB + C, (1.3a) 
«CC» = 2aT. (1.3b) 

This final theory is a stochastic Langevin equation, known as Bodeker's effective theory 
[3] or "the small frequency limit of Ampere's Law in a conductor" [4]. The parameter 
a is the "color conductivity." Theory 3 is only valid on spatial momentum scales small 
compared to the hard gluon damping rate 7 and frequencies small compared to momenta, 
cu< fc«7 = 0[s 2 Tln(<T 1 )]. 

Our goal is to calculate, by successive matching of these effective theories from short to 
large distance scales, the parameter a of Theory 3. 

The spatial scale at which physics becomes non-perturbative is k ~ g 2 T. So, by at least 
a logarithm, the interfaces (m, 7) of the successive effective theories are associated with 
perturbative physics, thus making a perturbative matching calculation possible. It will be 
useful to keep in mind a simple result from the analysis of static properties of hot gauge 
theories, which is that the parameter which controls the loop expansion is g 2 T/k, where 
k is the momentum scale of interest. So, in particular, the loop expansion for physics at 
the interface k ~ 7 between Theory 2 and Theory 3 is an expansion in inverse logarithms 

Our notation is the same as that of Ref. [1] and is summarized in Table I. W represents 
the adjoint color distribution of hard particles; SC is a linearized collision operator; £ and 
£ are Gaussian thermal noise; m is the leading-order Debye mass; and 7 is the hard gluon 
damping rate. Note in particular that (• • ■) denotes averages over the direction v of hard 
particle velocities, 1 whereas ((• • ■)) denotes averages over Gaussian noise. Also note that 
we use 5C to represent both a v-space integral operator 5C and the corresponding kernel 
5C(v, v'), which is simply a function. The formulas given earlier for the noise covariances 
((££)) and ((CO) should be understood as shorthand for 



Readers of ref. [5] should beware that our use of the notation (• • •) is completely different from theirs. 
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/=(l,v);va spatial unit vector. 

A = A(x, t), the spatial non-Abelian gauge field. 

W = W(x, v, t), the adjoint color distribution of hard excitations. 

£ = £(x, £) and £ = £(x, v, t) are Gaussian white noise. 

((• ■ ■)) denotes averaging over noise. 

(•••} = (•• -) v denotes averaging over the direction v. 

5 S2 (y— v') is a 5-function on the two-sphere normalized so that (5 S2 (v— v')) v / = 1. 
5CW = (5C(v, v') W^(v')) v /, the linearized collision operator applied to VT. 
" denotes either an operator on the space of functions of v, or a spatial unit vector. 
{ l m \ . . . \ Vm >) = J d n v Y^( v ) • • • Y lm (v) = 4tt(^(v) • • • Y* m (v)) v . 
(/| ■ ■ • |/') denotes the same for cases where the answer is proportional to S mtTn r. 
Pq = 1 00) <00 1 , the projection operator onto v-independent functions. 
7x = (1|5(7|1), the / = 1 eigenvalue of the linearized collision operator. 
Ca is the adjoint Casimir of the gauge group [N for SU(iV)]. 
d = 3 — e with e — > 0, the number of spatial dimensions. 
/' = /' d d q ^ r _ /"+ 00 rfgo 

iq ~ 7 (27T) d ' ig ~ 7-oo 27r' 

( — h ++) spacetime metric signature. 

A = q ,0 /q r , the ratio of frequency to spatial momenta. 

~ denotes equality at leading-log order. 

D = V + gA a T a , the gauge covariant derivative. 

T\ c = f abc , anti-Hermitian adjoint-representation generators. 

"In" in an order of magnitude estimate [e.g. 0(g 2 ln)] means ln(l/g). 

TABLE I. Summary of notation. 



2T 



((r(v,x,tKV,x',t')» = ^C(v,v'H a6 5 (3) (x-x') 



m 



<<C a (x,tR ? V,0>> = 2aT^,^5( 3 )(x-x') 5(t-t'), 



;i.4) 
;i.s) 



where i,j denote vector indices and a, b are adjoint color indices. 

The scale of the linearized collision operator 5C is set by 7. At leading log order, it is 
given by [3,4] 



5CW(v) = (5C(v,v')W(v')) v ,, 



;i.6a) 



5C(v, v') w 7 



,A2 



^(v-v')-- ,V ' V ' 



(1.6b) 



The symbol ~ denotes equalities valid only to leading-log order. 

We will consider the theories 1 and 2 to be ultraviolet (UV) regulated by dimensional 
regularization in d — 3— e dimensions with gauge coupling fi e ^ 2 g. Theory 3 is UV finite and 
does not require such regularization. 

In the remainder of this introduction, we review the path integral formulation of effective 
theories 2 and 3, and summarize general properties of the collision operator 5C which we 
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will need. In section II we perform the matching of Theory 3 to Theory 2, which is the 
most novel part of our calculation. The NLLO conductivity a is calculated in terms of 
the collision operator 5C of Theory 2. In section III, we determine the information we 
need about 8C by matching Theory 2 to Theory I. In the process, we explicitly calculate 
the hard thermal gauge boson damping rate in the presence of an infrared (IR) regulator 
(specifically dimensional regularization). Finally, in section IV, we put everything together 
and discuss the result. We then also summarize the differences between this work and an 
earlier discussion of the NLLO color conductivity by Blaizot and Iancu [5] . 



A. Review of Path Integral Formulation 

The original derivation of the Langevin Eq. (1.3) by Bodeker [3], and subsequent discus- 
sions [4], were performed in A = gauge, where the equation reads 

dA. 

-a— = BxB + C. (1.7) 

However, A = gauge is a sick gauge for doing perturbation theory, and is consequently 
an inappropriate choice for our present purposes. It is therefore useful to reformulate (1.7) 
as a path-integral, so that we can use standard Faddeev-Popov methods to choose a more 
convenient gauge. 2 Eq. (1.7) is a Langevin equation, and it is well known how to reformulate 
such equations as path integrals. 3 Specifically, Eq. (1.7) becomes 

Z = j [£>A(x, *)] exp {- J dt d 3 x Lj , (1.8a) 



with 

2 



l= 1 



AaT 



dt 



( j- 1 5 (3) (0)trD 2 , (1.8b) 



where D 2 means D-D. We will use dimensional regularization throughout our analysis, in 
which case one may ignore the trD 2 Jacobian term because 5^(0) = 0. Eq. (1.8b) is still 
in A = gauge, but we can now trivially generalize to a gauge-invariant formulation: 

Z = J [£>A)(x, *)] [£>A(x, t)} exp (- J dt d 3 x Lj , (1.9a) 



2 The Langevin equation (1.3) may also be shown to be correct and unambiguous in general flow gauges 
of the form Aq = R[A], where R[A] depends on A(x, t) only instantaneously and so does not involve time 
derivatives of A. See Ref. [6] for a discussion of flow gauges, and Ref. [7] for a proof that the equation (1.3) 
may be applied in any gauge of this class. There are subtleties, however, in directly interpreting the Langevin 
equation (1.3) in other gauge choices, such as Landau gauge. Using a flow gauge of the form A n = AV • A 
(which is discussed further in appendix A), it should be feasible to reproduce all the analysis of this paper 
directly from the Langevin equation. However, we found it most straightforward, both conceptually and 
computationally, to use the path integral formulations presented here. 

3 For a review, see chapters 4, 16, and 17, and in particular Eqs. (17.15) and (17.16), of Ref. [8]. See also 
Ref. [2]. 
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aE + DxB| 2 . (1.9b) 



AaT 



This can be checked by using the Faddeev-Popov procedure to return to A = gauge. But 
now we can use the Faddeev-Popov procedure on (1.9) to fix other gauges as well. Coulomb 
gauge, for instance, corresponds to 

Zcoioumb = J [VA ] [DA] [Vfj] [Drj\ 5(V • A) exp (- J dt d 3 x L Coloumb ) , (1.10a) 



L 



1 



Coloumb 



<xE + D x B\ 2 + f]V ■ ~Dr] , (1.10b) 



where fj and i] are anti-commuting Faddeev-Popov ghosts. 

Theory 2 can also be described by a path integral. We shall find it convenient to express 
the theory entirely in terms of the gauge fields by eliminating W using the equations of 
motion. The resulting path integral formulation is discussed in detail in Ref. [2], and here 
we simply quote the gauge-invariant result analogous to (1.9): 4 

L = -L[-ct(D)E + D x B] T a(D)- 1 [-a(D)E + D x b] +Li[A], (1.11) 

where <r(D) is now a matrix in vector-index space. It is also an operator in x space (and 
color) and is given by 

c^-(D) = m 2 Hm (v t [v • D + SC + AP ]"S) = m 2 {{viGvj) - fad) (G)' 1 (Gvj)) , (1.12) 

where G is the W-field propagator arising from (1.2a): 

G= [v-D + 5C}-\ (1.13) 

G is an operator in both x and v space (as well as color). In the first form of (1.12), P 
denotes the v-space projection operator that projects out functions that are independent of 
v. In the notation introduced below, Pq = |00)(00|. 

The term Li[A] in (1.11) is complicated and is discussed in Ref. [2]. It is the analog 
of the (^ 3 )(0) trD 2 term in (1.8b) but is spatially non-local and does not trivially vanish in 
dimensional regularization. Fortunately, however, the size of £i[A] is such that it will be 
irrelevant to our calculation of the NLLO conductivity. All we need to know about it for 
the present discussion is that £i[A] is independent of A , and that it is suppressed by the 
loop expansion parameter compared to other terms in (1.11). [Specifically, its contribution 
to the A propagator is suppressed by one factor of the loop expansion parameter.] As we 
will see later, this will be enough to argue that this term does not affect our calculation of 
the conductivity at NLLO. 



The path integral corresponding to this Lagrangian has time discretization ambiguities. These should be 
resolved by a time-symmetric prescription. That is, writing dA/dt = [A(t+e)— A(i)]/e, and then interpreting 
A's without time derivatives to mean the average [A(t+e)+A(t)]/2. 
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The trace of E ■ D x B is the space- time total derivative -V • (E a x B a ) - d t (±B a B a ). 
So, ignoring boundary terms, the cross-term may be dropped in the action for (l.ll): 5 

S = Jdtd 3 xL = Jdtd 3 x^ [Ea(D)E+(D x B)cr(D)- 1 (D x B)] + #%[A]j. (1.14) 

We should perhaps clarify our mixture of v-space operator notation and (• • •) notation 
for averaging over v. Since SC, Pq, and G are operators in v-space, they do not commute 
with Vi and Vj in (1.12). So, for instance, (vi 5C Vi) ^ (5C) even though ViVi = |v| 2 = 1. This 
particular example is made clear by rewriting (viSCvi) = (vi 5C(v, v') ^) V v'- Our notation 
works just like bra-ket notation in quantum mechanics if one rewrites (• • ■) as (00 1 ■ • • 1 00) 
with 1 00) representing the constant function Y 00 (v) = l/V^r in v-space. (It will be necessary 
later to consider other \lm) spherical harmonics as well.) 1 00) is only a v-space entity and 
does not specify anything about x or color dependence. So, for instance, (D) = (00|D|00) 
is D and not zero. 



B. General properties of SC 

We'll now discuss some useful properties of 5C that follow from general considerations 
not restricted to the leading log approximation (1.6b). Collisions packaged in 5C are, as far 
as effective Theory 2 is concerned, local in space. All of the x dependence of collision proba- 
bilities comes in the distribution functions W, and so the operator 5C itself is independent of 
x — it is simply an operator in v-space. Rotation invariance of the theory therefore implies 
v-space rotation invariance of the operator SC, which in turn implies that SC is diagonal in 
the space of |Zm)'s. [That is, its eigenfunctions are the spherical harmonics l^ m (v).] It also 
implies that the corresponding diagonal matrix elements (lm\5C\lm) depend only on I and 
not on m. We shall therefore write them more compactly as (l\5C\l). This /-decomposition 
of SC will be crucial to our later analysis, and the case / = 1 will be of particular interest. 6 

It's useful to demonstrate this notation by returning to the leading-log approximation 
(1.6b). The second term of (1.6b) vanishes when applied to any v-parity-odd function, and 
so (l\SC\l) ~ 7 for odd I. In particular, at leading-log order, 

(1|6(7|1)«7. (1.15) 

The relation is more complicated beyond leading-log order, and we will define a separate 
symbol 7 

71 MWl). (1.16) 



5 This reasoning depends on using a symmetric time discretization to define (1.11). Sec, for example, 
Appendix B of Rcf. [9]. Time is to be regarded as running from —00 to +00, and the choice of temporal 
boundary condition is irrelevant for our purposes. 

6 For a discussion of using the space of spherical harmonics \lm) as a basis for numerical simulations, see 
Rcf. [10]. 

7 In d dimensions, 71 means the eigenvalue of SC in the vector representation of SO(d). 
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In section II, we will find that the only pieces of SC that we need to calculate the NLLO 
conductivity are the leading-log formula (1.15), plus the NLLO result for 71. The eigenvalue 
71 will appear in our analysis through the handy formulas 

<JC'v)=7iv), (v5(7 = 7i(v. (1.17) 

Here, and henceforth, we will often use a bare ) or ( as convenient shorthand for 1 00) or (00 1, 
respectively. The relations (1.17) follow because v) has 1 = 1, and so is a superposition of 
|lra)'s. The value of 71 will be calculated in section III. 

Rotation invariance also implies that SC is symmetric in v-space, since 5C(v, v') can only 
depend on v • v'. Furthermore, like everything else in Eq. (1.2a), SC is real. As discussed in 
Ref. [2], this means that the conductivity operator <7 U (D) of (1.12) is real and symmetric in 
x/co lor /vector space. 

A very important property of SC is that it annihilates functions which are independent 
of v: 

SC) =0, (1.18) 

or equivalently 

(SC = 0, or (0|<$C7|0) = 0. (1.19) 

We will see this explicitly when we discuss SC in section III, but it is true for quite general 
reasons, as pointed out in Refs. [3,11]. One way to understand this is to observe that current 
conservation in effective theory 1 requires = = D^iv^W). Theory 2 is a subsequent 
effective theory for uj <C k, meaning that time derivatives have been neglected compared to 
spatial derivatives. In this limit, current conservation becomes simply (v-DW) = 0. Taking 
the v-expectation value of (1.2a) gives (v-DW) = (SC), and so (0015(7100) = 0. Since SC is 
diagonal in \lm) space, it follows that SC has a zero- mode, SC) = 0. A differently packaged 
but related explanation of this property may be found in sec. Ill C of Ref. [11]. 

Finally, the effect of a collision term in a Boltzmann equation is to cause the decay over 
time of correlations. The sign of SC in (1.2a) is such that decay corresponds to positive SC, 
which would be more obvious if we hadn't dropped the d t W time derivative term in going 
from (1.1a) of Theory 1 to (1.2a) of Theory 2. Based on this, one should expect that all the 
eigenvalues (l\SC\l) of SC are non-negative. We shall verify this explicitly in section II D 
using the leading-log formula (1.6b). The zero eigenvalue (1.18) corresponds to the fact that 
charge is conserved and does not decay. 

II. MATCHING THEORY 2 TO THEORY 3 

To match effective theories, one must identify gauge-invariant observables which are 
calculable in both theories. We shall therefore spend some time discussing a gauge-invariant 
observable, involving Wilson loops, that can be used to determine the parameter a of Theory 
3. We shall then find that, in practice, the matching problem can be conveniently simplified 
to the matching of the Coulomb-gauge self-energy n o of A in the limit of zero frequency 
and small momentum. 
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A. Wilson loops 



An example of a gauge-invariant observable that depends on the conductivity is a (real- 
time) Wilson loop 



that extends in the time direction, where V indicates path ordering of the exponential. Here 
= A ^ T a , and our convention is that the generators T a are ant i- Hermit ian. To see the 
dependence on conductivity explicitly, it is convenient to focus on rectangular Wilson loops, 
such as depicted in Fig. 1, where one set of edges is in the time direction and the other set 
is purely spatial. It will also be convenient to focus on rectangles whose temporal extent 
is very large compared to their spatial extent R. 



To get a feel for these Wilson loops, let's look at their value in our final effective theory, 
Theory 3, at first order in perturbation theory. There will be various perimeter contributions 
such as those of Fig. 2a. As we shall see, these are UV divergent and should in principle 
be regulated. 8 But they don't depend on the separation R, and so we can ignore them if 
we just focus on the independence of the Wilson loop expectation. We will similarly ignore 
contributions that don't depend on t^, such as Fig. 2b. R dependence is generated by 
propagators which connect different edges, such as in Fig. 2c. If we pick a reasonable gauge 
for doing perturbation theory, then, in the large time (too) limit for the Wilson loop, we can 
neglect diagrams such as Fig. 2d which attach to the far-past or far-future ends of the loop. 
In this case, the large-time Wilson loop, to lowest order, is determined just by Fig. 2c. 



8 As mentioned before, the dynamics of Theory 3 is UV finite. But the definition of the Wilson loop operator 
itself requires UV regularization. 




(2.1) 



R 



FIG. 1. A time-like Wilson loop rectangle. 



1. Relation to a in Theory 3 
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(a) (b) (c) (d) 

FIG. 2. Examples of first-order contribution to the expectation of a large-time Wilson rectangle. 

The primary example of an ■unreasonable gauge is A = gauge, which is the gauge in 
which Bodeker's effective theory (Theory 3) was originally formulated. There is no Fig. 2c 
at all in A = gauge, because the time-going Wilson lines only couple to A . But A = 
gauge is a sick gauge for perturbation theory in the first place. 9 We will instead work in 
Coulomb gauge. The calculation of this section is repeated in Appendix A in more general 
"flow" gauges, which interpolate smoothly between Coulomb gauge and A = gauge. 

The perturbative expansion of the action (1.10b) describing Bodeker's effective theory 
in Coulomb gauge is 



Scoioumb = Jdt d 3 x a\ VA | 2 + - I (ad t - V 2 ) A| 2 + 0(A 3 ) + (ghosts) 

One may read off the A propagator (which is instantaneous in time): 

IT ^ 



A a vAAAAA, A b = 



ak 2 



(2.2) 



(2.3) 



Fig. 2c then gives a contribution to the Wilson loop of 

ti(T a T b ) 



dn 1 SW = g 2 T[ ] jf °° dt dt' ((A a (t, 0)A b (t>, R] 



tr(l) 

= -g 2 C A t c 



2T 
k ak 2 



e ikR = 



2aC A T 
aR 



(2.4) 



where d R = tr(l) is the dimension of the representation associated with the Wilson loop. 
Because the A propagator is instantaneous in time, there are no crossed graphs at higher 
order and this contribution exponentiates in the usual way to give 



InW = 



2aC A T 
aR 



too + (higher order) 



(2.5) 



9 That's because the free action in A$ = gauge has infinitely many zero modes associated with time- 
independent gauge transformations — zero modes which are not properly treated in a perturbative expansion 
and which manifest as spurious non-intcgrablc singularities in propagators. 
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in the large-time limit, up to terms independent of R or t^. In perturbation theory, at least, 
one sees that a large time-like Wilson loop provides a gauge-invariant quantity from which 
one may extract a. 

One can automatically remove the perimeter terms independent of either R or by 
taking ratios of Wilson loops: 



In 



W(R 1 ,t 1 )W(R 2 ,t 2 ) 
W(R 1 ,t 2 )W(R 2 ,h) 



2aC A T f 1 1 
a \R[ ~ R~ 2 



(h - h) + (higher order). (2.6) 



This ratio is free of UV divergences in Theory 3. 

A warning is in order concerning the physical interpretation of the result (2.5). The result 
should be trustworthy whenever perturbation theory is reliable, which means whenever 
R <C l/g 2 T. Some readers may automatically associate Wilson loops with the behavior 
exp[—V(R)t 00 ], with V(R) interpreted as the potential energy (or free energy) associated 
with two static (i.e., infinitely massive) test charges separated by distance R. However, this 
interpretation only applies to Wilson loops (or Polyakov lines) in Euclidean time and does 
not apply to the case at hand of real-time loops at finite temperature. In particular, (2.5) 
should not be interpreted as a 1/R potential between static test charges for R <C 1/ g 2 T, 
which would be inconsistent with Debye screening for R 3> mr l = 0(1/ gT). In fact, we 
know no simple physical interpretation of real-time Wilson loops at finite temperature. 

Though we have only discussed a perturbative analysis of Wilson loops, that is good 
enough for matching Theory 2 to Theory 3, whose physics only differs at scales k > 7 = 
0(g 2 T In), where the physics is still perturbative (by a logarithm). Here and henceforth, in 
order of magnitude estimates the abbreviation "In" is shorthand for ln(l/g). As discussed in 
our companion paper [1], and in earlier works by Braaten and Nieto [12], matching may be 
performed by formally computing the same quantity perturbatively in both theories, in the 
presence of some common infrared regulator. We shall consider Wilson loops with R^> I/7, 
so that they are firmly in the region of validity of both theories, and we shall regulate the 
infrared behavior using dimensional regularization. For simplicity, we'll focus on the formal 
limit R — > 00 in the context of our IR-regulated perturbative calculation. (The order of 
limits is important: the — > 00 limit is to be taken first, so as not to invalidate the 
previous discussion.) 



2. Relation to a in Theory 2 

As far as the A propagator is concerned, the perturbative expansion of the action (1.14) 
is much like that of Theory 3 except that the color conductivity a becomes momentum 
dependent. Specifically, the Coulomb gauge propagator for A is now 

IT 

K vAAAAA, A\ = 5 ab , (2.7) 

k o"l (k) 

where 

^W^M^k)^-, (2.8) 
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and af) = o"ij| A=0 is given by (1.12) and (1.13) with the covariant derivative D replaced by 
ik: 



a^(k) = m 2 lim (^[v • ik + 5C + AP ]~ \ 

'(viGoMvj) - (^G (k))(Go(k))- 1 (G (k)^) , (2.9) 



m 2 



Here, 



G (k) = [v • V + 5C]- 1 = [v • ik + (2.10) 

The first-order contribution to the Wilson loop, analogous to (2.4), is then 

r 2T 

dn l = -g 2 C Atoo / - e* kR (2.11) 

One should expect that the large i? behavior is dominated by the small k behavior of the 
integrand, and this is indeed the case, up to corrections suppressed by powers of g. (See 
appendix B for an explicit argument.) The result is then 

(2 - 12) 

for R — > oo. More specifically, this limit is R ^> 7 _1 , since 5C provides the only scale in the 
definition of d"^ )(k) and the scale of 5C is 7. 
From (2.9), we have 

^ 0) (0) = m 2 lim ( Vi \5C + KP^v,) (2.13) 

As noted in section I, Vj) has 1 = 1 — that is, it is a superposition of |lm)'s. Recalling that 
5C is diagonal in the space of |/m)'s, as is P = 1 00) (00 1 , we obtain 



4°'(0) = ^, m d »f(0) = g, (2.14) 



in d spatial dimensions, where 71 = (1|5C|1) is the / = 1 eigenvalue of 5C. [We've left the 
spatial dimension d arbitrary in (2.14) because the generalization away from d=3 will be 
needed later when we dimensionally regularize.] Comparison of the Wilson loop (2.12) in 
Theory 2 and (2.4) in Theory 3 then gives the leading-order result for matching the two 
theories: 



2 9 
m m 

_ _ ... r~s^ 



d'-fi d^y ' ^ 

where the last leading-log equality uses (1.15). This is precisely the leading-log result for 
Theory 3 originally derived by Bodeker. This result could have been very quickly derived 
from the starting point of (1.6b) without all this discussion of Wilson loops. This approach 
based on Wilson loops does, however, provide a conceptually clear framework for discussing 
sub-leading corrections. 
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3. Next order in the loop expansion 



As mentioned earlier, perturbation theory at a scale k is controlled by the loop expansion 
parameter g 2 T/k. A matching calculation between Theory 2 and Theory 3 is a calculation 
of physics at the interface k ~ 7 below which Theory 3 is valid, and so the loop expansion 
for matching calculations will be an expansion in g 2 T/ r y ~ [lnf^l/g)] -1 . To go to next-to- 
leading-log order in the determination of a, we must therefore go to the next order in the 
loop expansion for the Wilson loops. 

A nice simplification occurs in Coulomb gauge. Consider the quadratic pieces of the 
Eff(D)E term in the action (1.14) for Theory 2: 

Ect(D)E = -A V • a(V)VA - 2A V • a(V) d t A - A a(V) d 2 A + 0(A 3 ). (2.16) 

Rotation invariance implies that Vj[<7(V)]jj must be proportional to Vj (since [o"(V)]y can 
only involve terms proportional to either Sij or VjVj). Therefore, the second term in (2.16) 
which connects Aq and A actually involves V • A and hence vanishes in Coulomb gauge. 
Consequently, with this gauge choice the propagator does not mix A and A. 

Now notice that the entire action (1.14) is quadratic in A — there are no Aq or higher 
terms. Taken together, this means that in Coulomb gauge there is no diagram such as Fig. 3 
contributing to the Wilson loop. Instead, the next-to-leading-order diagrams all have the 
form of self-energy corrections to the propagator in our leading-order diagram Fig. 2c. 



FIG. 3. A next-order correction to Fig. 2c that does not appear in Coulomb gauge. 

We shall henceforth visually distinguish A and A propagators when drawing Feynman 
diagrams, representing A propagators as dashed lines and A propagators as wavy lines. An 
expansion of the Theory 2 action (1.14) in powers of A will give interaction vertices of the 
forms shown in Fig. 4. In Coulomb gauge, the one-loop corrections to the Wilson loop are 
then given by Fig. 5. (We have left out diagrams involving tadpoles, as these vanish by CP 
symmetry.) This diagrammatic result applies to Theory 3 as well as to Theory 2. 




FIG. 4. The Coulomb-gauge interaction vertices of action (1.14) of Theory 2. Dashed lines represent 
Aq, wavy lines A, and solid lines the gauge-fixing ghost n. n is any positive integer. 
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(a) 



(b) 



o 



(c) 



FIG. 5. The first-order and next-order contributions to the Wilson loop in Coulomb gauge. 

In Coulomb gauge, the calculation of the Wilson loop at this order now reduces to the 
evaluation of the one-loop self-energy Hoo(u;,k) of Aq. In fact, for the large Wilson loops 
discussed earlier, all that will be important is the u=0, k — > behavior of the self-energy. 
Absorbing n 00 into the A propagator (2.7) and repeating the argument that led to (2.12), 
we have 



-(o) 



(k) 



(k) 



IT 

4\k) + — U oo (0,k), 



A: 2 



and 



d^ 1 5W 



2aC A T 



tr. 



(2.17) 



(2.18) 



in Theory 2. Corresponding expressions hold for Theory 3, with a^\k) replaced by a. 
Matching the two theories then gives 

^2 ( rOT 1 



m 

— h hm - 

(Z71 k^O 



2T 

-^n 00 (0,A;) 



(Theory 2) 



n 00 (o, k) 



(Theory 3) . 



a 



+ 0( J? ). (2.19) 



Here and henceforth, all results are in Coulomb gauge unless explicitly stated otherwise. 

It must be emphasized that, because this is a matching calculation, the k — > limits 
of the individual Hoo/A; 2 terms above are to be understood as taken in the presence of an 
infrared regulator. Suppose, for the sake of discussion, that Hoo(A;)/A; 2 were given by 

(2.20) 



k' 



P p 2 |k + p| 



where p was the spatial part of some loop momentum. The integral is perfectly finite and 
so does not appear to require any IR regularization. Without IR regularization, the result 
of (2.20) goes like k~ l by dimensional analysis, and the k — > limit is not well defined. Now 
imagine a simple momentum cut-off M on small loop momenta. The small k limit of (2.20) 
then behaves as M~ l rather than k~ l , and there is no problem with the limit. In general, 
in the presence of an IR regulator, we can formally expand integrands in k, so that, for 
example, 

in any IR regularization scheme. In dimensional regularization the result of this particular 
example is especially simple: (2.21) is zero. 
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B. Calculating n 00 (0,fc^O) in Theory 2 

1. Perturbative expansion of the action 

We will now be explicit about the perturbative expansion of the action (1.14) for Theory 
2. Since we are working in Coulomb gauge, there is an additional ghost piece of fj'V ■ Drj in 
the action but this is irrelevant since no ghosts enter any of the diagrams in Fig. 5 that we 
need to calculate. In fact, we will be parsimonious in our discussion of the expansion and 
only explicitly keep track of those terms which are relevant to the specific diagrams of Fig. 
5, namely the quadratic terms plus the A AA, A A A, and A A AA interactions. 

Before expanding the action (1.14), it is helpful to first make some simplifications. In 
particular, write E = DA — A and expand E<rE as 

j d d x EaE = J d d x [-A DaD A + 2A Da A + AaA] . (2.22) 

From the expression (1.12) for a, we have 

Bcr = m 2 ((v • D GV) - (v • D G) (G) _1 (GV)) . (2.23) 

Now comes a trick we will use repeatedly. Since 5C annihilates v-independent states, as 
discussed in the sec. IB, we have {5C = and can rewrite 

(v • D = ((v • D + 5C) = (G- 1 . (2.24) 

So 

(v-DG=(, and GvD) = ). (2.25) 
Eq. (2.23) then simplifies to 

Dcr = -m 2 (G')- 1 (G'v), (2.26) 

and 

DaD = -m 2 (G)-\ (2.27) 

We may now rewrite (2.22) as 

J d d x EaE = j d d x ^m 2 A (G)' 1 A - 2m 2 A (G)- 1 (Gv) ■ A + AaA] . (2.28) 

The perturbative expansion is now generated by formally expanding the expression (1.13) 
forG: 

G = [v • D + 5C}- 1 = G - gG Q w ■ AG + g 2 G Q v ■ AG v • AG , (2.29) 

with G given by (2.10). (We will not explicitly show the factors of /W 2 that accompany g 
in dimensional regularization.) Then 
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(Gy 1 = (Goy 1 + 5 (Go)- 1 (G'oV • AGoXGq)- 1 

+ ^(Go) -1 ^ • AGo)(G )" 1 (G'ov • AG ){G )- 1 
- /(G )- 1 (G v • AG v • AGo)(G )- 1 

+ 0(# 3 A 3 ) . (2.30) 

Now focus for a moment on the expansion of the A (G)~ 1 (Gv) ■ A term in (2.28), and 
consider in particular the leading-order contribution to (Gv) ■ A. By rotation invariance, 
(GoVi) must have a factor of Vj, since V is the only vector quantity appearing in Go- 
Combining this with 

(Gov) • V = (G (v • V + 6C)) = 1 (2.31) 

then yields 10 



(G v) = (vG ) = ^ • (2.32) 



Since we are in Coulomb gauge, we thus have (G v) • A = 0. So the leading term of the 
expansion vanishes, and therefore 

AoiG^iGv) ■ A = -gMGo^iGoV ■ AG v) • A + O(g 2 A A 3 ). (2.33) 

Putting everything together with the other terms in the action (1.14), and keeping track 
only of the terms that are needed for the diagrams of Fig. 5, yields S = Sf ree + S int , with 

Sirco = J dtd d x^m 2 A {G )- l A + Aa(V) A + AV^V)]" 1 V 2 A + (ghosts)}, (2.34a) 
5 int = |^^^|m%((G)- 1 -(Go)- 1 )A + 2^m%(Go)-XG v-AG v)-A 

+ (not needed)), (2.34b) 



where (G) -1 is expanded as shown in (2.30). The perturbative expansion of the Li[A] term 
of the action (1.14) falls into the "not needed" category since the diagrams of Fig. 5 do not 
contain any interaction vertices involving only A and not A , and any correction to the A 
propagator induced by L\ [A] will be of sub-leading order in logarithms. 

2. Propagators 

The perturbative a(V) appearing in Sf ree can be simplified a bit. In momentum space, 
we earlier called it a^(k), given by (2.9): 

(°) _ ™2//„, A „, \ _ / A 



= m'UviGoVj) - (^G )(G )- i (G ^) , (2.35) 



10 This formula does not easily generalize to (Gv) because of the non-commutativity of the covariant 
derivatives contained in G. 
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Now note that the first term is transverse, because 

ViiviGoVj) = ((V • v + 8C) G oVj ) = ( Vj ) = (2.36) 
(and similarly {v i G Vj)'S7 'j = 0). So we may rewrite 



m 2 ( Vi Go(k)^> = 4 0) (^)^(k), (2.37) 

where we introduce (perturbative) transverse and longitudinal projection operators 

P^'(k) = S ij - k l k\ (2.38) 
Pl j (k) = (2.39) 

may then be expressed as 

4o) (fc) = m2 k|o(kw_ (2 40) 

a — 1 

Similarly the second term in (2.35) is purely longitudinal, by (2.32). So 

(7<°>(k) = 4 0) (fc)Pr(k) + af (fc)ii(k), (2.41) 

with 



a? ) (fc) = ^-(G (k)>- 1 . (2.42) 



m 



2 



A; 2 

From (2.14) for the low-momentum limit a^°\0) we find 



2 



4 0) (0)=4 0) (0) = ^. (2-43) 

Because we are in Coulomb gauge, the longitudinal sector does not contribute to 
Aff(V) A or A[ct(V)] _1 A. So we can replace a by in the free action (2.34a). The 
resulting propagator for A is 



At v/VWW A) = S tJ (iv,k)5 ab ee — — 2 Pj?(k) 5 



The relations between these propagators and the retarded, equilibrium, or other types of 
propagators are discussed in Appendix C. The propagator for A Q is the same as (2.7) 
(although we are now representing this propagator by a dashed line), 

2T 2T - 

A o = 7^7777^ = 7^< G o(k)>^. (2.45) 

k 2 al (k) m 
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3. Transposition 

It will be convenient to be able to rewrite interaction terms (2.34b) as their transposes. 
Under transposition in x/color space, D T = — D, and so 

(G) T = ([v • D + 5C]- 1 ) T = ([-v ■ D + 5C]- 1 ) = (G), (2.46) 

where the last equality follows from taking v — > — v in the the v-average (•••)• The terms on 
the right-hand side of the perturbative expansion of (2.30) must also be symmetric, and this 
can be explicitly verified by recalling that the A's there are really color matrices A = A a T a , 
which are anti-symmetric because T h ac = f abc . So, for instance, 

(Gov • AG ) T = -(G^v • AGq) = (G v • AG ) (2.47) 

where the last equality again follows by v-parity and where we have implicitly used the fact 
that v and 8C are symmetric in v-space. 

Now let us check the transposition of the Aq^Gq)^ 1 (Gqv ■ AGqv) ■ A interaction in (2.34b). 
Being a little more explicit about color indices than previously, and placing an under-tilde 
on the A which is to be interpreted as a color matrix, the interaction is 

^(G'S b )" 1 (G'o c v- A cd G$ e v) ■ A e . (2.48) 

Its transpose is 

-A • (vGjv • AG r )(G r )- 1 A = —A • (vG v • AG )(G )A . (2.49) 
In summary, this interaction term can be written in either of two ways: 

A (G )- 1 (G'ov ■ AG v) • A = -A • (vG v • AG )(G )A . (2.50) 



4- Analysis of diagrams 



Instead of presenting general Feynman rules for all the various vertices in the effective 
theory, and applying these rules to the diagrams of Fig. 5, we will simply write the expres- 
sions for the loop diagrams directly by treating S int as a perturbation in the path integral, 
expanding the exponential exp(— <Si nt ), and explicitly taking all possible Wick contractions 
of the fields. For the case at hand, this is far more convenient. The loops (a-c) in Fig. 5 
represent the expressions 



-n, 



(a) 
00 



-n, 



(b) 
00 



g 2 m 2 " 
' AT 



gm 

It" 



(G )" 1 ((G'ov • AG )(G )- 1 (G v • AG ) - (G v • AG v • AG >) (Go)" 1 , 

(2.51a) 



(Gor^Gov ■ AGo>(G r%l <g> ^(Go)- 1 (G v ■ AG )(G )- 1 



(2.51b) 



tt(c) 



gm 
^XT 



(G y l f(G v • AG v • A) ® (-v • AG v ■ AG ) 



+ (G v • AG v • A) ® (-v • AG v • AG )1 (Gq)- 1 . 



(2.51c) 
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The ® notation above denotes where strings of color index contractions end. Specifically, 
(• • •) <S> (• • •) represents the color matrix (• • -) a (- • -) h , where a and b are adjoint color indices. 11 
Using the A propagator (2.45), one can rewrite (2.51b) as 



-n 



(b) 
00 



-2 



g 2 m T 
AT 



(Go)- 1 ^ • AG )(G )- 1 (G' v • AG )(G }- 1 



It neatly cancels the first term in (2.51a) to leave 



n 



(a+b) 
00 



rr(a) , n (b) 
ii 00 + ii 00 



g 2 m 2 
2T 



(Go)- 1 (G v • AG v • AG )(G )- 1 . 



(2.52) 



(2.53) 



The next step is to use the A propagator (2.44) and put everything into frequency and 
momentum space 



12 



n 



(a+b) 
00 



ng(o,fc) 



(0,k) 



ab 



1 1 

g z m z 
2T~ 

2 4 

g m 



ab 4T 2 



P0,P 



(Go^v^ G (p) vjT c db &(k)> 5 y (-po, k-p) , (2.54) 

P0,P 

(£ (kKT a c d G (p)u,) vl S n(~Po, k-p) %{po, p) 

(v 3 G (p) v^ b G (k)) - (^Go(k-p) v^ b g (k))] , (2.55) 



x 



where we find it convenient to introduce 



0o(k) = 



Go(k) 



(2.56) 

(Go(k)> 

Remember that we are interested only in the case of zero external frequency k . The loop 
frequency integrals are easy to do using the explicit form (2.44) of S^-, and give 

Ti?(q) 



/ Sij(-po,q) 



P 2 Sik(po,qi) Sji(-po,q2) 



2T 2 PJ r k (q 1 )Pi l (q 2 ) 



9l^T 0) (92) + 92 2 4 UJ (?l) 



2=(0), 



Simplifying the color factors, we then have 

C A g 2 m 2 



n 



(a+b) 
00 



ng(o,fc) 



p^(k- P ) 



(0, k) = / p (^o(k) Vi G (P) ^ &(k)> ^ 



(2.57) 
(2.58) 

(2.59) 



G A £ 2 m 4 



/ (^o(k)^Go(p) 



(k-p) if (p) 



k-p| 2 4 0) (p)+p 2 ^(|k- P |) 



X 



(wjGo(p)^o(k)) + (v T Go{k-p)vjQo(X)) 



(2.60) 



n So, for example, (G v- AG v- A)® (v AG v- AG ) represents (Gg c v- A cd G^ e v- A e ) (v-A^G^v- A^G^). 

12 The relative minus sign in (2.55) arises from the time derivatives in (2.51c), which give (— ipo)(—ipo) f° r 
the first term, and (— ipo)(+ipo) for the second. 



5. Thek^O limit 



We now want to extract the k — > limit of H 00 (k) through Oik 2 ), as this is required 
for the Wilson loop matching formula (2.19). This section is somewhat tedious, and some 
readers may wish to skip to the result (2.81). 

To classify the various terms in the k expansion, it will help to rewrite the total self- 
energy given by (2.59) and (2.60) as 



n 00 (0,A;) = ^(£o(k)d(k)£o(k)>, 



(2.61) 



with 



6(k)=m 2 \v i G (p)v 



|k— p| ; 



-m 2 Vi G Q (p) vj 



k p A p 



|k-p| 2 o- p +p 2 o- k _ p 
where in this section we use the abbreviations 

- ~(°) 



(ujGo(p) % + (v T G (k-p) ■ 



P« = P?(p), 



[2.62) 



(2.63) 



to keep formulas more compact. The operator O is symmetric in v-space. This is manifest 
for all the terms except the one involving Go(k— p), and it is easily seen for that term by 
making the change of integration variable p — > k— p. 

It is useful to start with the small k expansion of Q (k)). This expansion is derived in 
appendix D and is 



&(k)> = 



(Go(k)> V 7i J 1 



[2M) 



We now turn to the operator 0(k), and first consider 0(k)). From the definitions (2.62) of 
0(k) and (2.37) of a p = cr^\p), we have 



6(k)> = m 2 / v t G (p)v 



gVg K + q k _ p ) 



and in particular 



|k-p| 2 |k-p| 2 (T p +p 2 o- k _ p 



O(0)) = 0. 



^2.65) 



(2.66) 



This means that 0(k)) is O(k), and it is the reason that we only need the expansion of 
Qo(k) thru 0(k) and not 0(k 2 ). Specifically, using the symmetry of Oik) and Qo(k) under 
v-space transposition, we can organize the small k expansion of (2.61) as 



^oo(k) = 



C A9 2 



(£>(k)) -2 7l - 1 (vikO(k)) + 7l - 2 (vikC>(0)v-ik) +0(k 3 



(2.67) 
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First term: Let's begin with (0(h)). Starting with (2.65), one can again use the defi- 
nition (2.37) of Op = (Jt\p) to find 



(0(k)) 



pij pV n 
r p Mc-p°P 



(Cp + °"k-p) 



|k- P | 



/ pij pV 
/ ^P^k-p 
Jp 



(p 2 ~ |k- P | 2 ) 
|k-p| 2 



|k-p| 2 (Tp +p 2 (Tk-p 
OpOk-p 

|k-p| 2 o- p + p 2 o- k _ F 



(2.68) 



Now symmetrize the integrand with respect to the change of integration variable p — > k— p: 



(0(k)) = I P^Pl p 
= (d-1) f 

Jp 



(p 2 - 


k-p 


) o-p^k-p 


2p 2 \ 


k-p| 


2 (|k-p 


2 o- p + p 2 a k _ p ) 



(k-p)Vp 



+ 0(k 3 ) 



(d-1) 
d 



a 



k 2 



p T 



(2.69) 



Second term: For the r ) 1 x (v -ikO(k)) term of (2.67), we need the 0(h) piece of 0(k)}. 
Its extraction can be simplified by rewriting (2.65) as 



6(k))=m 2 [ ViGo^Vj 

Jp 



|k- P | 



Mc P Mc-p^p 

|k-p| 2 -p 2 )o k - P 



p«fc pjk 

|k-p| 2 (|k- P | 2 o- p + p 2 o- k _ p ) k p p t 



m 2 J v i G (p)v j ) \ 



P 



m 2 f Vi G (p) Uj) 4 fcV - 2p^p • k + <^p • k + 0(k 2 ) . 



P 



(2.70) 



Therefore, 



(v • ik 0(k)) = m 2 / (viVi G (p) Uj) — fcV - 2pV P • k + 5 ij p ■ k + 0(k 3 ) . (2.71) 



Because of rotation invariance, (v • «kC?(k)) depends only on the magnitude k of k, so 
nothing is harmed by averaging the integrand in (2.71) over the direction k, giving 



(vikO(k)) 



m 2 k 2 



1 



d Jpp A 



(Go(p)v-ip) -2(v-pv-pG (p)v-ip) + (v-ipuj G (p) u, 

(2.72) 



The factors of Go(p) can be eliminated from the first and second terms by using the trick 
(2.24), which in the present context is 



Go(p)v.ip) = G (p) (v-ip + SC)) 



(2.73) 
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Using the relations (1.17), a similar manipulation shows that 

(v-ip^Go(p) = (uiV-ipGo(p) = (vi(y ■ ip + 8C - 7l ) G (p) = (^[1 - 71 G (p)] , (2.74) 
and so the final term of (2.72) becomes 



7i 



(v • ipv i G (p)v i ) = 1 - (d-1) — <7 P . 



Putting everything together, 



7 r 1 (v-ikc»(k)) 



(d-l) 



A: 



pP 



(2a - (J p ) , 



(2.75) 



(2.76) 



where <t = m 2 /d'~f 1 is the value of a p at p = 0, taken from (2.43). 

Third term: Finally, we pursue the 7f 2 (v • ik O(0) v ■ ik) term in the expansion (2.67). 
We are again free to average over the direction k, giving 



(vikO(O)vik) = -—(v k O(0)v k ). 



(2.77) 



Now use the definition (2.62) of 6(k) to find: 



(v k 6(0) v k ) 



m 

~~2 



p V 



(vk Vi G (p)vjV k )Pg 



-m 2 (v k Vi G (p)vj 



pit pjj 
P P 

2ffn 



(vjG Q (p)viV k ) + (v T G (-p)vjV k ) \. (2.78) 



The second term vanishes by the following symmetry argument. First, use v — > — v and the 
definition (2.10) of Go(k) to rewrite 



{vjG {p)viV k ) + (v- t G (-p)vjV k ) = (vjG (p)v T v k ) - (1 <-> j). 



(2.79) 



By rotation invariance, (vjGq(p) VjV k ) can only depend on the direction as pipjpfc, piSj k , 
pj5j k , and p k 5jj. Every one of these possibilities either has a p, or a pj, which will annihilate 
against the transverse projections PpP£* in (2.78), or else has a 5 l \ which vanishes by the 
anti-symmetry of (2.79) in ij. In summary, (2.77) and (2.78) become simply 



/ 2 2 

(v • ik O(0) v • ik) = -- [ ?L (v k Vi G (p) V j v k )P*. 



[2.80) 



Simplifying this expression will require us to be a little more systematic about the manipula- 
tions we have been using and is the subject of a later section. We end this one by combining 
the results (2.69), (2.76), and (2.80) for the individual terms appearing in (2.67): 



n o(fc) 



c A9 2 e 



^-^ (Aa - 3<r p ) + ( Vk Vi G (p) Vj v k )P^ 



dp 4 



7ir 



+ 0(k 3 ). (2.81) 
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6. Extracting UV and IR divergences 



As we shall see, the integral (2.81) giving n 00 is both infrared and ultraviolet divergent 
in three spatial dimensions. We are using dimensional regularization, but it will simplify 
our discussion of what to do with (v k V{ Go(p) v j v k) if we can instead work directly in three 
dimensions. Therefore, we will now isolate the pieces of the integrand responsible for the IR 
and UV divergences, so that we can subtract them and evaluate the remainder as a finite 
integral in d— 3. Specifically, we will rewrite (2.81) as 



C A g 2 k 

n oo (*0 = tt~ 



0(k 3 ), (2.82) 



f /reg(p) + / /m(p) + / /UV(P) 

P J p J p 



where 



/reg(p) = ^TT (4a - 3a p ) + (v k v l G (p)v J v k )P^ - f^\p) - /( uv )(p). (2.83) 
dp 71P 

/ir and /uv will be chosen to (a) make the / reg integral finite, and so evaluable directly in 
three dimensions, and (b) make the /ir and /uv integrals analytically tractable in dimen- 
sional regularization. 

IR behavior: In appendix D, we show that Go(p) has the small p expansion 



d 



G = J (7i - v • ip) P (7i - v • ip) + 0(p°). (2.84) 

7iP 



We can thus expand the (vvGqvv) term in (2.81) as 



-^{v k Vi G {p) VjV k )P l J = ^p- (v k Vi (71 - v • ip) )((7i - iv • p) VjV k ) P' J + 0(p 2 ) 

= ^p +0 (p-). (2.85) 

As discussed earlier, the p — > limit of <t p is o"o = m 2 /{d r )i). In Appendix D, we show that 
the small p corrections to a(p) are 0(p 2 ). Putting everything together, we may then choose 



/lR (p) = 2^p. (2.86) 



In dimensional regularization, the integral of /ir vanishes: 

/ Zir(p) = 0. (2.87) 
J P 

UV behavior: For p — > oo, we can treat 8C as a perturbation to v • ip, giving G — > 
(v • ?p) _1 , except that we will need a prescription for integrating over the pole v • p = 
in angular integrals. The prescription is obtained by recalling from section I B that 5C is a 
non-negative operator. So 

G (p) = P-P. A- + 7r5(v • p) , (2.88) 

v • ip + e v • ip 



22 



where e is a positive infinitesimal and P.P. denotes principal part. (Higher-order corrections 
to this formula are discussed in Appendix E.) This limit then gives 



m 



(d-1) 



(v i G (p)v, l 



m 



(d-1) 



7T(5(v-p)) 



m 2 S d _ 2 7T 



and 



(v k Vi G (p) VjV k )I™ -> 7r(ui 5(v • p)uj)i^ = tt(5(v • p)) 



<Sd-2 



(2.89) 



(2.90) 



where 



Sd-i = 



2ir d / 2 
V{d/2) 



(2.91) 



is the surface area of a (d— l)-sphere (e.g., Si = 2n and S2 = 47r). The UV piece of the 
integrand (2.81) is therefore 



, / \ 5d_2 7T(7o 1 

JuvlPJ — ► -5 o 



(2.92) 



and comes only from the (vvG vv) term. We don't want /uv to mess up our IR subtraction, 
so we will cut it off in the infrared by choosing 



/uv(p) 



Sd-2 na 



1 



Sd-i7i P(P 2 + M 2 )' 
where M is arbitrary. In dimensional regularization, the integral of /uv is 



(2.93) 



Axv(p) 



T (e/2) (j / n 



2(4vr) 1 - e / 2 71 \M/ 4^ 
where /x is the MS scale defined by 



In 



JL 

M 



+ 0(e) 



(2.94) 



M = AM 



/4tt 

e 7B 



(2.95) 



Putting everything together, Eq. (2.82) for Iloo(fc) becomes 



n 00 (A;) 



Ca«o"o ^2 



271 



+ In 



ii 

M 



+ 



4-7T7 



O"0 



+ 0(A; 3 ) + 0(e) 



i / / reg (p,M) 

where we may now set d—3 in 

Wp. «) = £ (!-» - 2 -p) + ^5 <<* »< Go(p) «i «*> ^ - 271 p(p2 + M2) 

This formula for n 00 does not depend on the choice of M. 



TTCTo 



1 



(2.96a) 



(2.96b) 
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7. Reducing (v h ■■■v im G (p) v h ■ ■ ■ v jm ) 



In the previous section, we encountered several v-averages of G'o(p) flanked by various 
factors of v. Any finite combination of v's can be rewritten as a superposition of spherical 
harmonics |/m)'s, and so we can recast the problem of simplifying general expressions of the 
form (fjj • • -v im Cro(p) Vj 1 ■ ■ -Vj m ) to the simplification of (l'm'\G (p)\lm) . 

Choose the z-axis in the direction of p. Then 



(l'm'\G (p)\lm) = (l'm'\(ipv z + 5Cy l \lm). 



(2.97) 



Now recall that 8C is diagonal in / and to, and note that v z can change / but does not change 
the azimuthal quantum number to. So 



(l'm'\G (p)\lm) = (l'm\(ipv z + 5C)~ l \lm) 5 mm >. 

We can derive a recursion relation in I by writing 

5w = (l'm\(ipv z + 5C)Go\lm) 

= ip(l'm\v z \l"m) (l"m\G \lm) + (l'\5C\l') (l'm\G \lm) , 

with an implied sum over I". Since v z can only change I by ±1, this gives 



(2.98) 



(2.99) 



8 W = ip(l'm\v z \(l'+l)m)((l'+l)m\G \lm) + ip(l'm\v z \(l'-l)m)((l'-l)m\G \lm) 

+ (l'\6C\l'){l'm\G \lm) . (2.100) 

This defines a recursion relation which allows one to rewrite matrix elements with higher 
/' in terms of those with lower I'. This recursion will end when /' becomes as low as it can 
be consistent with m — that is, at V = \m\. A similar recursion can be constructed for /, 
and by use of these recursions, all matrix elements {l f m\Go\lm) can be rewritten in terms of 
{\m\m Gq \m\m). In fact, the case m < is related to the case m > by v-parity: 



(\m\,m Go(p) \m\,m) = (\m\,—m Gq(— p) 



ml 



m) = (\m\, — m Go(p) \ m \, —m) , (2.101) 



where the final equality follows from v-parity and the [v • ip + SC] 1 structure of Gq(p). 

We now turn to the specific problem of rewriting the {v^Vi G'o(p) v j v k) term of our ex- 
pression (2.81) for n o in terms of the (mm\Go\mm) . A product of two v's is a combination 
of 1—2 and 1=0, so we will be able to rewrite the expectation in terms of (22|Go|22) and 
(Go). The advantage of this rewriting is that later analysis of how to evaluate expecta- 
tions involving G will be simpler and more natural for (mm\G \mm) than for (vkViG VjVk) 
directly. 

It is easiest to work backwards from the explicit form for F 2 ,2 which gives 



|22) = -1 



32tt 



V X + Wy) 



(2.102) 



From this, one may easily check that 



24 



(22|G (p)|22> = ^ (v k v t G (p)v jVl } (2P^P kl - ifpfj , 



8 



(2.103) 



for any operator G (p), remembering that we have chosen p to be the z direction for the 
purpose of defining \lm). 13 From (2.103), we have 



(vkViGoVjVk)!™ = {v k v i G v j v l )P; j (P kl + p k p l ) 



= — (22|G |22) + (vuvAvjVi) (^^' + J?P*p') • 



15 

The second term on the right-hand side can be simplified by expanding P^> 
and by repeated use of the relation (2.74). The result is 

(vkViGoVjvjP* = ^ (22|G |22) + \ (G 0/ , 
The analogous derivation for (ll|(?o|ll) gives 14 



7i (\ 2 Op 
6 3 ctq 



(2.104) 
(2.105) 



;il|G |ll) = I (viGovJP? = \ (v t G v t ) = 



(2.106) 



To obtain the conductivity at NLLO, it is adequate to use leading-log approximations 
to the propagators in our one- loop calculation of Hoo- That is, we only need to evaluate the 
integrand / reg in (2.96) at leading-log order. At this order, if one factors out the scale 7 of 
5C, functions like (22|G f o(p)|22) may be re-expressed as purely numerical functions of the 
single dimensionless variable p/7. Specifically, at leading-log order, 



m\ m 



G (p) |m|m)«7 1 T,\ m \(p/l), 



where 



S m (p) 



mm 



1 



mm 



ipv z + 5c 

and 5c is the leading-log result (1.6b) for 5C '/j: 

4 (v • v') 2 



(2.107) 



(2.108) 



8c(v,v') = 5 S2 {v- v') 



(2.109) 



13 The natural generalization to d dimensions is 

(22|G (p)|22) = { J^_ 2) (v k Vi G (p)v jVl ) [(d-1) P^P? - P^Pg] , 

where the relative coefficient of the two P P P P terms is chosen so that contraction with 5ik or Sji gives zero 
(so as to exclude the I = combinations of vv), and the overall normalization has been chosen so that 
replacing G (p) by 1 gives (22|22) = 1. 

14 Thc natural generalization to d dimensions is (11|G |11) = — — (u, G a Vj) P^ = da p /m 2 . 
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We can now combine this with (2.96) and (2.105) to obtain 



C A aa 2 
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+ In 




x 



1 + (in" 1 )] +0(k 3 ) +0(e), (2.110) 



with 



/(!/) = - / dp 
7T JO 



2 ^ 



_3_ A _8_/ 
(?) + 3p 2 V 



15 ^ 



7rp 



2(p 2 + z/ 2 ) 



(2.111) 



We have used (2.106) to rewrite a p (at the order under consideration) as <7 £i(p/7). We 
have ceased to distinguish between 7x and 7 in n o since we are ignoring corrections to n 00 
suppressed by additional powers of inverse logs [see (1.15)]. The terms in (2.111) have been 
arranged so that each term is individually IR safe. We will discuss how to evaluate I{y) 
numerically in section II D. 



C. Matching to Theory 3 



We've now got n o in Theory 2, but we still need n o in Theory 3, so that we can 
use the matching condition (2.19) to determine the parameter a of Theory 3 at NLLO. 
Fortunately, the one-loop calculation in Theory 3 is trivial in dimensional regularization: 
lim fc ^ [^ -2 n o(A;)] = 0. The reason is simple dimensional analysis. Rescale the variables of 
the path integral (1.9) for Theory 3 to t = a'H, A = T~ 1/2 A, and A = aT- 1/2 A . Here 
we will for once be explicit about the factors of pf. The path action can then be rewritten 
as 



dtd d x 



-E + D x B 



(2.112) 



where D u = d v + gpfl 2 T x l 2 A v . In this form, the parameters of the theory appear only in 
the combination gpf^T 1 / 2 . At one loop, the self-energy n 00 = er _1 Tn o of A must be 
proportional to g 2 p, e T, which has mass dimension 1 + e. But lim fe ^ [k~ 2 fl 00 (A;)] has mass 
dimension zero, and there are no other dimensionful parameters in the problem that can 
make up the discrepancy in mass dimension! Consistency then forces limfc^ [k -2 Iloo(^)] 
= in dimensional regularization. Such simplicity is the standard virtue of dimensional 
regularization for matching calculations [12,13]. 

Our matching condition (2.19) and our Theory 2 result (2.110) then give the color con- 
ductivity a at NLLO: 



a 



m 
<*7i 



7 



+ o(hr 



(2.113) 



where, for the sake of definiteness, we have fixed M — 7. As one can see, the only information 
we will need about 8C at NLLO is the value of 71. 

We now turn to methods for evaluating the functions S m (p) and so evaluating the nu- 
merical constant 1(1) from (2.111). 
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D. Evaluation of S Tn ,(p) « 7 (mm|Go(p)|mm) 



The dimensionless functions X m (p) were defined in (2.108) as {mm\(iv z p + 5c)' 



\mm) 



where 5c is the leading-log result for SC/j. As discussed in section II B 7, the operators v z 
and 5C both preserve m, 5C preserves I as well, and v z changes I by ±1. For fixed m, the 
operator iv z p + 5c may therefore be considered as a tri-diagonal matrix in the \lm) basis 
where I = m, m+1, m+2, ... : 



iv z p 



5c 



( 



Cm 
m P 



c m +i 



*6 (m) 



to m+lP 
C m +2 

ih {m) n 

W m+2P 



W m+2P 



Cm+3 



\ 



(2.114) 



where 



q = (/|5c|/>, &[ m) = (Zm|u z |(Z+l)m). 



(2.115) 



£ m (p) corresponds to the upper- left element of the inverse of the matrix (2.114). Inverting 
tri-diagonal matrices is particularly simple, and there is a continued-fraction formula for this 
element: 15 



S m (p) 



1 



Cm "i 



(2.116) 



Cm+l + 



(6 {m) 



m+l 



pY 



C m +2 + 



Note that X m (p) is even in p. 

All we need now are explicit formulas for the coefficients q and b\ . Eq. (2.116) may 
then be used for numerical evaluation of X m (p). 16 The b\ m ^ are given by 



b\ m) = 



{I + l) 2 - m 2 
^ 4(/ + 1) 2 -1 



(2.117) 



15 This easily follows from iterating the formula for the inverse of an (N+l) x (7V+1) matrix in terms 

of the inverse of its lower-right N x N block: ^ ^ ^ ^ = ^ ad~ 1 c (\+ad^ >1 cb T ) d^ 1 ) ' w ^ a ^ 
(a — b T d _1 c) _1 . Here, a is a scalar, b and c are TV-component (column) vectors, and d is an TV x N matrix. 

16 From Eqs. (2.117) and (2.119), one may read off that &j m ' — > ^ and c; — > 1 as I — ► 00. Consequently, the 

tail of the continued fraction, X = lim^oo c m+ k + (&|"+fc/o) 2 / '(■ • •), satisfies JV = 1 + (p/2) 2 X~ 1 , implying 

that JV = 5 + 5\/l + p 2 - The continued fraction (2.116) may be evaluated quite accurately by replacing its 
tail at large but finite I by this value. 
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The q may be evaluated from the expression (2.109) for 6c(v,v') as 



q = (5c(v,v')PKv-v')>, 



l- 2 -f +1 dz z2p ' (z) 



7T J-l v' I - -2 

where -Pz(z) are Legendre polynomials. The integral vanishes if / is odd and gives 

1 2 



(2.118) 



17 



C 2 r. 



1 - 2 



(2n)! 
2 2n {n\f 



1 + 



(2n) 2 + 2n - 2/ ' 



C2n+1 — 1 • 



(2.119) 



The q are all non-negative, as was claimed in section IB. Note that Co vanishes, as it must. 

The procedure for numerical evaluation of S m (p) is to compute the continued-fraction 
formula (2.116) with some upper cut-off / max on /, and then repeat the calculation, doubling 
Z max each time until the answer converges. This procedure becomes inefficient for very large 
p, however, because it then requires rather large / max for good convergence. For very large p, 
it is more convenient to use asymptotic formula for £ m (p), which are derived and presented 
in Appendix E. 

The final result of numerical evaluation of the integral (2.111) that defines 1(f), using 
numerical evaluation of the functions £ m (p) as described above, is 



1(1) = 2.8380' 



(2.120) 



III. MATCHING THEORY 1 TO THEORY 2 

Our next task is to determine the operator SC of Theory 2, appearing in (1.2a). Specif- 
ically, we want 71 = (115(711) to leading order in g (and all orders in logs). We will follow 
the general matching strategy used previously. We will temporarily introduce an infrared 
cut-off, then compute the total effective collision operator 5C tot in both theories, formally 
expanded to leading order in perturbation theory, and then determine what bare collision 
operator 5C appearing in Eq. (1.2a) of Theory 2 is required for the results to match. We 
will again use dimensional regularization to regulate the infrared. 



A. SCtot in the underlying theory (Theory 1) 

There are now a variety of methods for computing the effective collision operator at 
leading order in the underlying short-distance theory [3, 15, 4, 16, 17]. 18 Previous authors 



17 We evaluated the integral using Eq. 2.17.2 of Ref. [14], and verified the result numerically. 

18 Another interesting analytic approach that gives 71 at leading log order is that of sections 2 and 3 of ref. 
[18]. It is not clear to us how to extend this approach beyond leading log order, and in particular how to 

obtain the terms of 5C(v,v') that are not proportional to 5(v — v'). 
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have only extracted the leading log piece of their result because the leading-order result 
is formally log divergent in the infrared. Having regulated the infrared, we shall instead 
extract the entire thing. 19 So, one may now follow, in d spatial dimensions, one's favorite 
method of the references just cited. The method we're most intimately familiar with is our 
own, so our discussion will most closely parallel the presentation in Ref. [4]. 

At leading order in g, 5C is generated by 2 <-> 2 collisions of hard particles, mediated by 
semi-hard (qo < q ^ m) t-channel gluon exchange, such as depicted in Fig. 6. One finds 20 

5C w{w) = (I l " M(v ' v '' q) |2 [w{v) " ^ (v ' )] ) v , ' (3 ' 1} 

where M. is the amplitude for a t-channel collision between hard particles with velocities v 
and v', mediated by a semi-hard gauge boson with momentum q. (This interpretation of M. 
reverts to a more fundamental picture than that of Theory 1, interpreting W as made up of 
individual hard particles. A derivation that is more directly in the framework of Theory 1 
may be found in Ref. [3].) 



FIG. 6. The dominant scattering process: t-channel gauge boson exchange. The solid lines represent 
any sort of hard particles, including the non-Abelian gauge bosons themselves. 

If the integrand of (3.1) is separated into two pieces, the coefficient of the first term, 
proportional to W(v), is the same expression one obtains in a leading order calculation of 
the hard gauge boson damping rate [19], 

^w(/,^ (v ' v ' q)1 l' (3 ' 2) 

The overall coefficient in front of / q 1 | 2 in (3.2) simply represents the results of group 
factors and the integration of the magnitude |p'| in the calculation of 7 based on Fig. 6. 
The dependence on exactly what species of hard particles are present is completely isolated 
in the value of the Debye screening mass. The formula (3.1) for 5CW(v) may equivalently 
be converted to a formula for the kernel 5C(v, v') itself, 



xr>f >\ C A m 2 T 
06 v, v = 



/ |M(v,v',q)| 2 \ «f*(v-v')- / |A4(v,v',q) 

Jq / v / Jq 



(3.3) 



19 Ref. [15] also discusses doing this but does not go so far as to extract an explicit result. 

20 Here, as well as in Eqs. (1.1) and (1.2) defining either effective theory, m should be understood as the 
d-dimensional Debye mass. 
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We have seen that, for the calculation of the NLLO conductivity, we will not need the 
full form of 5C, but will only need the matrix element 71 = (1|5C|1) = (v l 5Cv l ). We will 
begin with a general analysis of SC, however, and specialize to 71 later. 

If there were no screening effects to consider, the scattering amplitude M would be the 
classic Coulomb amplitude such that 



gY e 



f d 4 ~*Q 


6^ , 


1 (2tt) 4 -^ 





2tt6(Q ■v)27r6(Q-v'), (no screening) (3.4) 



where Q = (q°,q) and v = (l,v). The (^-functions are simply q°,q C T approximations of 
the constraints that the final-state hard particles be on shell. To account for screening of 
the exchanged semi-hard gluon, however, we must replace this by 



gY" 



(2tt) 4 - 



Pt U (Q) 



+ 



PT(Q) 



^g 2 + n T (Q) g 2 + n L (g) 



2tt5(Q ■ v) 2n5(Q ■ v') , 
(3.5) 



where Pp and Pl are the transverse and longitudinal projection operators, 

0, \x = or v = 0; 



T(g) = 



5 ij - 



gT- 



q l q 3 
q A 

w 
g 2 



\i = % and v = j, 



(3.6) 
(3.7) 



and where IIt and I1l are the transverse and longitudinal pieces of the standard, leading- 
order, hard thermal loop self-energy [20]. This self-energy depends on Q only through the 
ratio A = q°/q. 

We will review the explicit formula for 11(g) below, but it's worth first examining some 
qualitative features. The longitudinal sector is screened for small Q by Debye/plasmon 
effects. The transverse sector, however, is unscreened in the A — > limit, reflecting the fact 
that charged plasmas do not screen static magnetic fields. This lack of static screening is 
responsible for the logarithmic infrared sensitivity that generates the usual leading-log result 
for the color conductivity. The logarithmic divergence appears only in the purely transverse 
term of the squared amplitude (3.5). It will be convenient to isolate that divergence by 
rewriting (3.5) as 



(\M\ 2 - \M\ 



\M\ 2 m , 



(3.8) 



where |.M| 2 R is a small q° (small A) limiting form of \M\ 2 that we shall discuss in a mo- 
ment. Our strategy is to arrange that the integral of |.M| 2 R be simple enough to evaluate 
in dimensional regularization, whereas the first term in (3.8) will be completely finite and 
evaluable directly in d = 3 dimensions. 

To continue, switch integration variables from q° to A and, making use of the ^-functions, 
rewrite (3.5) as 
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/ \M\ 2 = 1vg^ 

/q 



, + 1 r 1 
l - 2 ' dX ' 



q q 



v • v' - A 2 



+ 



-1 + A 2 



g 2 (l-A 2 )+n T (A) g 2 (l-A 2 )+n L (A) 
x (J(A-q-v) 5(A-q- v') . 



(3.9) 



The angular integration over q is equivalent to replacing the pair of delta functions by their 
angular average. We will denote this average, in d- dimensions, by 



f d (X, v • v') = (5(X - q • v) 5(A-q.v')).. 
We will implement our split (3.8) by extracting the small A behavior of Eq. (3.9), 

r+l P 1 ,r . 2 

|2 



(3.10) 



\M\i R = 2ng' l fi 



, +i 

/ dX 



-i 



1 
q q 



v ■ v 



g 2 + n^ R (A) 



5(q-v) 5(q-v'), 



where is the limiting small A behavior of Ht, to be discussed explicitly below. 



(3.11) 



1. The IR piece 



In terms of 



f d (0, v • v') = <5(q • v) tf(q • v') 



we have 



4 2tt^/1 - (v • v') 2 ' 



/ |^| I 2 R = 27r(v.v') 2 /,(0,v.v')//i 2e / + ^A / 



(3.12) 



qg|g 2 + n^ R (A)r 



(3.13) 



We now need the form of n! ] R (A). The usual d = 3 result for the small frequency behavior of 
ITt is — ^ m 2 A. But we need the result in 3 — e dimensions. One could derive this directly by 
evaluating a one- loop thermal diagram in the fundamental quantum field theory, but let's 
instead derive it in the usual way from effective theory 1. Working at leading order in g, 
formally solving the W equation (1.1a), and plugging into the Maxwell equation (1.1b) gives 



d v F^ v = m 2 



v ■ d + £ 



(3.14) 



where the e is simply an infinitesimal prescription specifying retarded behavior. Comparing 
with d v F^ v + TV^A,, = 0, the components of the self-energy can then be read off, among 
which 



U ij (q) = -m 2 q° 
We want the transverse self-energy 



v v- 



-q° + v • q - is / v 



Tin 



d-l 



tr(P T n). 



(3.15) 



(3.16) 
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One can evaluate this for arbitrary frequency, 21 but here we're only interested in the small 
frequency limit nJj?\ Taking q° small, (3.15) and (3.16) yield 



n 



IR 



m 2 q° 1 1 — (v ■ q) 



d— 1 \ v • q — is 



m 2 q° 
' d-1 



(in5(v ■ q)), 



(d-l)S d . 



where the area S^-i of a (d— l)-sphere is given by (2.91), and where 

2 



K 



e + 0(e 2 ). 



-5(i,|)=l+[ln2- 

We can now do the momentum integral in (3.13) using 

S d -i vr/4 



m A = — - nm A, 
4 

(3.17) 



(3.18) 



/ 1 

Jqq\q 2 + ic\ 2 Jqg(g 4 - 



I (rf— 5)/2 



+ c 2 ) (27r) d cos(e7r/4) 



(3.19) 



to give 



/ \M\ 



2 

IR 



r+1 



rfA 



T^T — T~7IY ( v ' v ') 2 ^(°> v • v ') ^ r 

{27i) d 1 cos(e7r/4) J-i 

(v- v) /40,V'V)g /i (^-rem J 



—nm A 
4 



-l-(e/2) 



-l-(6/2) 



e(27r) d - 1 cos(e7r/4) 
2#V(v • v') 2 
n 2 m 2 ^Jl — (v • v') s 



+ In 



/i 



I l _ i 

~r 2 |E o 



(3.20) 



5. T/ie non-infrared piece 



Now turn to the remaining piece, / d-M] 2 — |-M|ir)> which may be evaluated directly in 



three dimensions. From (3.9) and (3.13), we have 



M\ 2 -\M\ 2 m )=2ng 



+1 d\[ - 

-l Jq q 



v • v' - A 2 



+ 



-1 + A 2 



g 2 (l-A 2 ) + n T (A) g 2 (l-A 2 ) + IT L (A) 



v • v 



A2 



|g 2 + n^ R (A)|' 



/ 3 (0,v.v') . 



The explicit form of / 3 is 
/ 3 (A,v-vO 



2tt 



1 - (v • v') 2 - 2A 2 (1 - v • v') 



-1/2 



/ 3 (A,v- v') 
(3.21) 

(3.22) 



(This is really multiplied by a step function which vanishes when the argument of the square 
root goes negative.) The basic 3-dimensional q integral required is 



The result is n T = ^ [A 2 + (1-A 2 ) 2 Fi(±, 1; f ; A~ 2 )] and n L = m 2 (l-A 2 ) [l - 2*1(5, ^ f 5 A ~ 2 )] ■ 
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1 In Z2 — In Z\ 

q g(g 2 + zi)(q 2 + z 2 ) A-n 2 (z 2 -zi) 



(3.23) 



where the cut of the logarithm is understood to run along the negative real axis. It is 
convenient to rewrite Hl and Ht as 



IIl,t(A) = m 2 (1— A 2 ) Pl,t(A). 



We then obtain 



/ q (^i 2 -i^iy = ^/> {/,(a,vv 



vv'-A 2 V 



1-A 2 
v • v' — 



1-A 2 



IA 



/ 3 (0,vv') (vV) 



/\2 



(3.24) 



argp T (A) argp L (A) 
Imp T (A) Imp L (A) 



A 2 \ Re lnp^(A)-lnp L (A) 



Pt(A)-Pl(A) 



(3.25) 



where argz = Im (In z) is to be understood to lie in the range [— ir, n]. 

Finally, we need explicit formulas for pt and in three dimensions. One can look up 
the formulas for II [20] or easily derive them from Eq. (3.15). In the case at hand, we are 
interested in space-like momenta Q, and the results are 



Pl (A) 
p T (A) 




(3.26a) 
(3.26b) 



By combining Eqs. (3.3), (3.8), (3.20), (3.22), (3.25), and (3.26), we now have a complete, 
if somewhat cumbersome and inelegant, integral formula for 8C(v, v') at leading order in 
g. Because of the remaining A integration, the functional dependence of SC on v • v' is not 
simple. Fortunately, we do not need the complete form of 5C(v, v') to calculate the NLLO 
conductivity, and we will now specialize to the calculation of the matrix element 71. 



3. Calculation of 71 and 7 
Using the general formulas (3.1) or (3.3) for 5C, we have 

7iH^CV> = %^(/ |M(v,v',q)| 2 (l-vv')) . (3.27) 

2#V Vq / vv / 

As it turns out, we can easily calculate 7, defined by (3.2), at the same time as 71. So we 
will, even though we don't actually need the NLLO value of 7 for our calculation of the 
NLLO conductivity. To this end, we define 

7 (r ' ) = %S?(/l^(v,v',q)| 2 (l-r ? v.v')) , (3.28) 

zg fj, Vq / vv / 
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where 77 = 1 yields 71 and 77 = yields 7. We now apply this to the pieces (3.20) and (3.25) 
of / \M.\ 2 using the three dimensional identities 



/ 3 (A,v-v')(l-r ? vv / ; 



S(\ — q • v) S(\ — q • v') (1 — 77V • v'' 



= J(i-»A 2 ), 



qvv 



and similarly 

((v ■ v' - A 2 ) / 3 (A, v • v') (1 - t?v • v')) vv , = -lv (1 " A 2 ) 2 , 

((v ■ v' - A 2 ) 2 / 3 (A, v ■ v') (1 - 77V • v'))^ = |(1 - A 2 ) 2 (1 - r]\ 2 ) , 

plus the d = 3 — e identity 

(v • v') 2 / d (0, v • v') (1 - 77V • v')) , = ((v • v') 2 5(q • v) *(q ■ v') 



(3.29) 

(3.30) 
(3.31) 



qvv' 



's, 



d-2 



Inserting Eq. (3.20) into (3.28) then gives 



7 £> = -C A aT 



C^aT 



2^-2 



d-1 \5 d _i 



1/771 



(3.32) 



(d— l)7r5' d _iC0s(e7r/4) ne \4p/ 
-I + lnQ +I 7E -21n2 + 0(e) 



7T \ e 2 



And inserting the non-infrared piece (3.25) into (3.28) now produces 

1 W- 1 $=C A aTaW, 
where the numerical constant is given by the one dimensional integral 



(3.33) 



(3.34) 



" lA|}' 



1 argp T (A) argp L (A) 

2 Impx(A) ImpL(A) 



+ 7?(l-A 2 )Re 



lnp T (A)-lnp L (A) 
p T (A)- PL (A) 



(3.35) 



It is useful for numerical evaluation to split pt and pl into their real and imaginary parts, 
Pt = Rt + i-^T, e ^ c - ; an d use 



arg p 
Im p 



cot~ 



(3.36) 



and 



Re 



Rt—Rl 



lnp^(A) - In p L (A) = 1 

p^(A) - p L (A) ~ 2 [(it> T -it> L ) 2 + (J t +/l) 2 ] "V^l + 



ln /i? 2 +/f 



|/t+/l| 



[(.Rt-^l) 2 + (/t+/l) : 



cot 



-if -^L 



cot 



-1 f -fr] 



(3.37) 
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where cot" 1 a; = | — Tan^x is defined to lie in the range [0, 7r], and we have made use of 
the fact that the signs of Jl and Jl + It are the same and are both opposite to It- One may 
also note that the integrand of (3.35) is even in A. 
Putting everything together, we find 



and 



71 = C A aT 



7 = Ca olT 




i In - + oi 

4 



— + ln — + Mn- + a 

e 2 4 



(3.38) 



(3.39) 



where we have now written the result in terms of the MS scale jx = \^bte 7E / 2 /i . Numerical 
evaluation of (3.35) for r\ = 1 and rj = gives 



and 



^ = 0.323833- 
a = 0.120782 • 



(3.40) 
(3.41) 



for T) — 1 and rj — 0, respectively. In fact, our numerical evaluation of the constant a shows 
that it precisely equals — |ln| to 12 significant digits. Surely this is an exact identity, 22 
so that the dimensionally regulated hard gauge boson damping rate, to next-to-leading-log 
order, is simply 



7 = C A aT 




(3.42) 



B. Matching to Theory 2 

In theory 2, the total effective collision term is, in principal, composed of two parts. 
First, there's the bare collision term that appears in (1.2a), which we will call <5Cba re here 
to be explicit, and which conceptually represents collisions due to the exchange of virtual 
gauge bosons that were integrated out in going from Theory 1 to Theory 2. Secondly, there 
is a dynamical contribution to SC, which we will call 5Cd yn , which arises from the exchange 
of those gauge bosons that have not yet been integrated out. However, as we explain 
below, the same nice property of dimensional regularization which simplified the Theory 3 
matching calculation in section II C works here as well: for the purposes of matching, 8C^ yn 
must vanish in dimensional regularization by dimensional analysis. Hence, we have simply 
$C tot = 5Cb are , and so the 71 we needed in order to match Theory 2 to Theory 3 can simply 
be taken directly from the Theory 1 result (3.38) for (l|<5C'tot|l)- 



22 We have not bothered to try proving this analytically, although we are confidant one could do so. Blaizot 
and Iancu [5] have shown the corresponding result when using a sharp momentum IR cutoff, instead of 
dimensional regularization. (See appendix B of [21].) 
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The dimensional argument can be made by rescaling to variables t 



m 



H, A = T-V2A, 



A = m^-VMo, 
Theory 2 become 



W = m 2 T- l l 2 W, and f = m 2 ^ 1 / 2 so that the equations (1.2) for 



v • D W - v • E = -5C bare W + l 
(W) = 0, 
D x B = (vW), 

((£>, V, i) V', t')» = 2 5C ba *c(v, v') (J(x-x') *(*-*'), 



(3.43a) 
(3.43b) 
(3.43c) 
(3.43d) 



where £>„ = 9 V + g^l 2 T x l 2 A v . There remains a dimensional quantity, 5C, other than the 
effective coupling g[fl 2 T x l 2 . But in matching theories 1 and 2, it is important that <5Cb are is 
to be formally treated as a perturbation. This is feasible because the matching is performed 
at the spatial momentum scale m = O(gT), whereas SC^, are is of the much more infrared 
scale 7 = 0(g 2 T In). As discussed earlier, matching can always be thought of as taking 
place in a large box, serving as an infrared regulator. The box should be chosen to be large 
compared to the distance scale of the matching (1/m), but may be small compared to more 
infrared scales where the physics becomes more complicated (e.g., I/7). In the presence of 
such an IR cut-off, one may then treat quantities that are soft relative to the matching scale 
(such as 5Cb arc , in the case at hand) as perturbations. Moreover, formally treating them as 
perturbations works for the purposes of matching calculations even if the infrared regulator 
is then taken to arbitrarily large distance scales [12,13]. 

In summary, then, perturbation theory in the effective coupling gpi € l 2 T x l 2 and in 5C can 
only give terms with the dimensions of {^g[fl 2 T x ^ (j*C^J for integer m and n. Because of 

the factor of //, none of these can consistently match the mass dimension 1 of SC tot unless 
m — 0. That means the only contribution to 5C tot is the tree-level tfCWe- 

Though the dynamical contributions to 8C formally vanish in Theory 2 in dimen- 
sional regularization, they do play an important conceptual role. We are now interpreting 
Eq. (3.38) as the result for the bare 71 in Theory 2. Theory 2 requires ultraviolet regular- 
ization: the 1/e in (3.38) is the counter-term for a UV divergence in Theory 2, and fi is 
the associated renormalization scale. In our matching calculation, however, the 1/e and the 
ln(m//x) actually arose in Eq. (3.20) from a formal infrared divergence of the calculation of 
the total 71 in Theory 1. The discrepancy of interpretation is resolved by realizing that the 
Theory 2 result that 71 tat — 7i,bare should really be thought of as 



[7l,tot] (Thcory2) - 7l,bare + C A aT 




+ 



IK 




(3.44) 



and so equating the total 71 in the two theories converts the the IR divergence into a 
UV divergence in (3.38). Eq. (3.44) is an example of the generic behavior in dimensional 
regularization of logarithmic divergences when there is no scale to cut them off in either the 
IR or UV, and is typified by the simple example 



d 3 ~ 



p 



OC /I 



dp 



p 



l+e 



dp 



p 



l+e 



dp 



A p 



l+e 



1 , A 

h In — 

e ji 



+ 



1 , A 
H In — 

e pi 



(3.45) 



where in the last equality dimensional regularization was used both for the IR contribution 
of the first term and the UV contribution of the second. 
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IV. FINAL RESULTS 



We now put together our NLLO result (2.113) for the color conductivity a in terms of 
71 and our result (3.38) for 71. The structure is clearest if we write an expansion for a 1 
(the "color resistivity" ) rather than a directly. One finds 



_ 1 3C\aT 
a = 



m 2 



ln (^)) +c+0( " rl) 



(4.1a) 



with 



C= |ln^ + ai + 7(1) = 3.0410- •• , (4.1b) 



and where the numerical constants 1(1) and a\ are defined by (2.111) and (3.35) [with i]=l] 
and given numerically by (2.120) and (3.40). Note that the 1/e divergences have canceled, as 
they must. Inside the logarithm of (4.1a), 7(71) is to be understood as simply the leading-log 
formula 

7 (//) « C A aT\n Q , (4.2) 

and \x should be chosen so that it is of order 7. One may easily verify that the fi dependence 
in the NLLO result (4.1a) only affects that answer at order [ln(m/7)] _1 ~ [ln(l/g)] _1 , which 
is beyond the order of this calculation. 

Eq. (4.1) is our final result. Although the result must be gauge independent, our deriva- 
tion has been restricted to a particular choice of gauge, namely Coulomb gauge. It would be 
comforting to have the calculation repeated in another gauge, perhaps the generalized flow 
gauges discussed in appendix A, but this we have not attempted to do. It would also be 
interesting if there were any way to express the fundamental functions S m (p) of this problem 
in terms of standard mathematical functions, but we have been unable to do so. 

Our result for the NLLO color conductivity may be compared against numerical simu- 
lation of the electroweak baryon number violation rate, as such simulations can in fact be 
used to measure the relative size of the NLLO correction to a. We discuss this comparison 
in ref. [1]. 

Finally, we should mention the differences between our calculation of the NLLO conduc- 
tivity and that outlined in earlier work by Blaizot and Iancu [5]. In our language, Blaizot 
and Iancu's discussion amounts to specifying how to calculate 71 (which they call 7 + 5). 
That is, it is equivalent to our discussion of matching Theory 1 and Theory 2, though they 
do not actually push through the calculation to get a final, regulated, numerical result like 
(3.38). They then quote a result for the conductivity of simply m 2 /(37i). However, this 
misses all the contributions that in our calculation came from integrating out physics at 
k ~ 7 when matching Theory 2 to Theory 3. Moreover, the result m 2 /(37i) is not a well- 
defined quantity at NLLO because of the infrared divergence in the calculation of 71. If 
one simply cuts off that divergence at k ~ 7, as suggested by Blaizot and Iancu, 71 is still 
sensitive at NLLO to one's convention in choosing whether to cut it off at k = 7, k = 27, 
or something similar. The calculation of the effects treated in our Theory 2 to Theory 3 
matching is crucial to obtain an answer that is independent at NLLO to the precise choice 
of cut-off scale (renormalization scale) /i. 
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APPENDIX A: WILSON LOOPS IN FLOW GAUGES 



A useful class of gauges for stochastic gauge theory are the flow gauges defined by the 
condition aA = — AV • A. This class of gauges interpolate between A = gauge (A — > 0) 
and Coulomb gauge (A — > oo) [6]. 23 In this appendix, we will illustrate the use of these 
gauges by explicitly checking the gauge-invariance of the derivation in section II A of the 
the first-order result for large-time rectangular Wilson loops. 

The gauge-fixed path integral, analogous to (1.10), is 



with 



4w = J [DAq] [DA] [Vfj] [Vrj] 5(aA -XV • A) exp (- J dt d 3 x L 
L flow = ^ {\-aE + D x B| 2 + fj(aD + AV • D)^} . 



(Ala) 



(Alb) 



Rewriting A in terms of V- A and working in momentum space, the perturbative expansion 
of the action gives 



Show = J dt 
The propagator is 



AaT 



(ad t - V 2 ) A + (l-A)VV • A| 2 + 0(A 3 ) + (ghosts) 



(A2) 



A a vAAAAA/ A b = 2aT5 ab 

1 3 



\iauj + k 2 \ 2 liauj + Xk 2 \ 2 



(A3) 



One can see that the gauge A = 1 is the stochastic gauge theory analog of Feynman gauge. 
This Feynman-like gauge was first proposed in this context by Zwanziger [22]. 

In any case, let's now turn to the Wilson loop diagrams of Fig. 5c. We need the propa- 
gator for A = (A/(t)V ■ A which, from (A3), is 



IT 

An ^VWW A b n = —5 



ab 



\ 2 k 2 



\iouj + Xk 2 



2 • 



(A4) 



Then, following (2.4), 



= ~9 2 C A t c 



2T 



\ 2 k 2 



c \iauj + AA; 2 1 



e ik-R = 



J ui=0 



aC A T 
aR 



(A5) 



in the large-time limit. As required for any physical quantity, all dependence on the gauge- 
fixing parameter A has disappeared. 



21! 



The flow gauges here correspond to the £ — > oo limit of the generalized flow gauges considered in Ref. [6] . 
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APPENDIX B: LARGE R BEHAVIOR OF WILSON LOOPS 



Consider the first-order correction to the real-time Wilson loop expectation (2.11), re- 
peated here for convenience, 

5W = -2g 2 T C A f «, ( /!* , (Bl) 
J * k 2 crl J (k) 

in the limit of large R. Specifically R must be large compared to the color-changing mean 
free path 7 -1 = O [{g 2 T In) -1 ]. Physically, one should simultaneously keep R sufficiently 
small that physics on the scale R is still perturbative, that is, R <C (g 2 T)~ l . Formally, 
when performing an IR regulated matching calculation, one may simply take R — > 00. But 
here, we will examine what happens if one is not quite so cavalier with R. For the sake of 
definiteness, we will consider R ~ 5(g 2 T)~ 1 , where l/lng^ 1 <C 5 <C 1 and the coupling g is 
arbitrarily weak. 

Rewrite the Fourier transform in (Bl) as 

m) = / k i? e * k ' R > ( B2 ) 

where s(k) = [cr^^A;)] -1 , and perform the angular integral: 

I poo fjh 

™ = *niL T* k) **"">■ (B3) 

The analytic continuation of s(k) is an even function of k, as can be verified explicitly from 
the formulas (2.42) and (2.116). So (B3) can be rewritten as 

^ - h ^-r r dk sJ v [ e " s - e ~ aR \ - ^ l dk ' kR rr (f) ■ < B4 > 

Here, P.P. denotes principal part, 

PP. (\) !(^' -.'!. (B5) 



k) 2 \k — ie k + ie 

Closing the contour of (B4) in the upper half plane then picks up the k = ie pole, as well as 
any contributions from singularities in the function s(k). The k = ie pole gives the expected 
result (2.12). Singularities in s(k) = [o'^\k)]~ 1 at some k in the upper half plane will give 
additional contributions suppressed by exp[— RIm(k)\. a^\k) is analytic and non- vanishing 
in a neighborhood of the real axis. The only scale associated with the k dependence of af\k) 
is the scale 7 appearing in 5C. So the only scale that can determine the imaginary parts of 
the singularities of s is 7. Therefore, the contribution to (B4) from singularities of s must 
be suppressed by exp[— 0{^R)\. For R ~ 5{g 2 T)^ 1 this is exp[— 5 0(lng~ r )] which, for weak 
coupling, behaves as some (positive) power g and vanishes faster than any power of 1/ In g~ x . 
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APPENDIX C: RELATIONS BETWEEN VARIOUS CORRELATORS 



Consider, for simplicity, an effective theory for a single real scalar field, governed by 
the classical Langevin equation a<p = (V 2 — m 2 ) — dV int /d(fi + ( and the noise covariance 
((CO) — 2<tT. Also for simplicity, take a to be constant. The corresponding action in a path 
integral formulation is then 24 



1 



S = — I dt d 3 x 
4aT 



cr0 - (V -m)(f>- 



(Cl) 



If SV/S(p is treated as a perturbation, then the unperturbed propagator arising from this 
action is, in Fourier space, 

((<f){uj, k) <f){uj', k')*» = -iG(u, k) (27r) 4 5(k-k') 5{uj-uj') . (C2) 

with 

2aT 

-iG(u,k) = - — ^ (C3) 

In the main text, the term "propagator" referred what we've called — %G above, which in 
coordinate space is ((<f>(t, x) 0(0, 0))). In the underlying quantum field theory, however, it is 
customary to define the propagator G as i times the expectation of a product of fields, as 
indicated above. Note from (C3) that G is purely imaginary. 

Because operators commute in a classical theory, there is no distinction between time- 
ordered and time-unordered propagators. So the above propagator can be regarded as the 
classical limit of the time-ordered propagator 25 G(t,x) = i(T{0(t,x) ^(0,0)}), or either of 
the Wightman functions G>(f,x) =i(^(t,x) </» f (0,0)) or G<(t,x) = i(0 f (O, 0) <j)(t, x)), of the 
underlying quantum theory. That is, 

G = G > = G < (classically). (C4) 

A classical effective theory, such as the Langevin theory (CI), is (at best) valid for low 
frequency, long wavelength dynamics. It is only appropriate for studying observables which 
are smeared over time and spatial scales large compared to (3K, and thus insensitive to short 
time (or short distance) quantum fluctuations. 

In the underlying quantum theory, the spectral density p(u, k) (defined as the Fourier 
transform of ([(j>(t, x), (f>'(0, 0)])) is related to the time-ordered propagator G by 

p(u, k) = 2 tanh (\(3fiw) Im G(u, k) . (C5) 



24 Formally, the action should also contain the Jacobian term — \ tr(— V 2 + m 2 + V^ t (</>)) . But this is 

proportional to 8^(0) for local potentials Mnt(0) and so vanishes in dimensional rcgularization. See, for 
example, sec. 17.2 of ref. [8]. 

25 In our specific example, (f> is real field, and so ${t, x) is no different from <j)(t,x). In writing general 
quantum expressions for propagators and spectral densities and so forth, we nonetheless find it useful for 
the sake of reference to include daggers where they would appear for a discussion of complex fields. 
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Although it is not essential to the discussion, we will keep track of factors of h in this 
appendix (and this appendix only). Taking the classical limit hou <C T, in which the real- 
time propagator G is purely imaginary, (C5) gives 

ip(w,k)=phuG(uj,k). (C6) 

The Fourier transform of the retarded propagator GrR,(i,x) = iO(t) ([(f>(t, x), <ft(0, 0)]) is 
related to p by the spectral representation 

G R(Wlk)= / + ~^^L, (07) 

J-oo 2-K UJ' — UJ — IE 

which, in the classical limit, becomes just 

G R (t,k)=i9(t)(3h^G(t,k). (C8) 

Similarly, the advanced propagator G/i(t,x) = —i9(—t) {[4>{t, x), 0^(0, 0)]) satisfies the same 
spectral representation (C7) except for changing — ie to +ie, and in the classical limit is 

G A (t, k) = -i9(-t) Ph ^ G(t, k) . (C9) 

Note that these propagators satisfy the standard relation G R — Ga — ip. 

In studying static equilibrium physics, one typically works in imaginary time r 
rather than real time t, where Euclidean correlation functions such as Ge(t, x) = 
(T{0(— it, x)0t(O, 0)}) are periodic with period (3h. In the classical limit, which is the 
long-distance limit Ax 3> (3h, there is no substantive difference between equal-time correla- 
tions and zero-frequency correlations in imaginary time, due to the negligible extent (3% of 
imaginary time. In this limit, the only discrepancy is a change in normalization of (3% from 
the Fourier transform, 

G E (r=0, k) = {f3hy l G E (u=0, k) . (C10) 

The limit of zero time separation is the same in imaginary or real time, so one must have 
G(t=0, k) = 2Ge(t=0, k), since the real-time time-ordered propagator is related to the Eu- 
clidean propagator by appropriate analytic continuation in t. But the analytic continuation 
in frequency of Ge from the Matsubara points v n = 2mn/ {(3%) back to real frequencies yields 
Gr or Ga, depending on whether one continues to the real frequency axis from above or be- 
low, respectively. Hence, the zero frequency retarded, advanced, and Euclidean propagators 
coincide, GE(^=0,k) = G R (c<j=0,k) = Ga(co—0, k). Putting everything together gives the 
following string of equalities in the classical limit: 

G E (r=0,k) = {Phy 1 G E (v=0,k) = (Ph)- 1 G R (uj=0,k) = (/^GA^Ck) 

/+°° duo ~ 
— G(u,k). (Cll) 

Note that this also implies 



42 



/ + 00 f + OO 

dt G R (t, k) = i(ph)- 1 / dt G A (t, k), (C12) 
-OO J —CO 

which can be seen directly from (C8) and (C9). 

It is important to understand that the small frequency limit of the real-time correlator 
G(oj, k) is not directly related to the zero frequency limit of the Euclidean correlator Ge(^, k). 
As an example, for our real-time propagator (C3), 

^=°- k > = (FT^r (C13) 

while 

r+co jT 

G(t=0, k) = j ^ dco G{u, k) = . (C14) 

In this simple scalar model, the chain of relations (Cll) imply that Ge{t=0,]<.) and 
(Ph)' 1 Ge(V=0, k) likewise equal T/(k 2 + m 2 ). Some additional specific results for this 
model are 

k ) = 7 — \2 2^2 ' ( C15 ) 

(auy + (k 2 + m 2 ) 2 

G R (^k) = , \ (C16) 
— icru) + k 2 + m 2 

G A (co,k)= , -A— - . (C17) 
zero; + &r + 

Finally, if one inserts the spectral density (C15) directly into the spectral representation, in 
the underlying quantum theory, of the Euclidean propagator 



M=r^, (cis) 

J-oo 27r a;' + zi/ n 



then one obtains 



^ k ' %H^^ - (C19) 

However, this result is not sensible (for non-zero Matsubara frequencies). As noted earlier, 
the classical Langevin theory which led to (C15) is only valid for frequencies small compared 
to (Ph)" 1 , whereas the result (C19) depends sensitively on the precise form of the high 
frequency behavior of the spectral density (C15). 



APPENDIX D: SMALL p EXPANSION OF G (p) 

In this appendix, we discuss how to expand Go(p) — [v ■ ip + SC] -1 for small momentum 
p. We cannot simply treat iv ■ p as a perturbation because 5C has a zero mode (1.18) and 
so is not invertible. The correct way to expand is to separate iv • p into zero-mode and 
non-zero mode pieces, 

v • ip = v • ip P + P v • ip + (1-P ) v • ip (1-P ) , (Dl) 
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and treat only the last term as a perturbation [11]. The first two terms are rank-one operators 
which will lift the zero mode of SC. This leads to the expansion 



Go= [G 00 1 + (l-Po)v-2p(l-Po) 

= Goo - G 00 (1-Po) v • ip (1-Po) Goo 

+ G 00 (1-P ) v • ip (1-P ) G 00 (1-Po) v • ip (1-P ) G 00 - 



(D2a) 



where 



Goo = 



v-ipP + P v-ip + <5G 



d 



(1-Po) 5C' 1 (1-Po) + — (71 - v • ip) P (71 - v • ip). 



7iP 



(D2b) 



To verify the last equality, note that SC v • ipPo = 71 v • ipPo by (1.17). Observe that Goo 
is 0(p~ 2 ). That might appear problematical for the expansion (D2a), which brings along a 
factor of Goo with every factor of v • ip, except that the inner factors of Goo always appear 
in the combination 



(1-Po) Goo (1-Po) = (1-Po) 



5C 1 — — v • ip P v • ip 

7i 



(l-^o), 



(D3) 



which is only 0(p°). 

The expansion (D2) yields 



Go 



d 



7iP 



(71 - v • ip) P (71 - v • ip) + 0(p°), 



and so 



and 



G > = ^(7i-v-ip)) + O(pP), 



(Go) = ^- + 0( P ° 



(D4) 



(D5) 



(D6) 



Putting the last two equations together gives the expansion (2.64) for Q ) cited in the main 
text. 

The expansion (D2) also gives 



m 

a *=d^l 



{viGoVi) = a + 0{p 2 ) , 



(D7) 



which is a result used in section (II B 6). 
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APPENDIX E: LARGE p EXPANSION OF E m (p) 

Using the large p expansion (2.88) of G (p), one finds 

f 7T 

^m(p) ~ (mm|P.P. h ir5( y v z p)\mm) = — (mm\ 5{v z ) \mm) 

ID % P P 

(2m+f) 2 tt (2m+l)!!(2m-l)H vr 
-72^)T [Pm(0)] Yp~ (2m)\ 27' (E1) 

where (— 1)!! = f. [The P.P. terms vanishes because the \mm) states are invariant under 

We can make a little progress analyzing the corrections to this leading term by attempting 
to treat 5c as a perturbation in [iv x p + e) -1 , 

S m (p) = (mm\(iv z p + e) _1 |mm) — (mm\(iv z p + e) -1 <5c (if z p + e) _1 |mm) + • • • (E2) 

Call the second term of the expansion 5S m . Taking the formula (2.109) for 5c, we can rewrite 
this term as 



5Yi m (p) = (mm\(iv z p + e) 2 \mm) 



-47T^ m (v) 

where 



P.P. — + vr5(^p) 
iv z p 



5c 2 (v, v') 



P.P. t^- + 7T(5(^p) 



F mm (v')) , (E3) 



vv' 



4 (v • v') 2 

5c 2 (v,v') = r — ■ (E4) 

^^l-(vv') 2 

By independently considering both v — > — v and v' — > — v' in the second term of (E3), one 
sees that only the 5-functions will contribute for even m and only the principal-part for odd 
m. 

As we shall see, the vv' average in the second term of (E3) is logarithmically divergent 
as simultaneously v z — > and v' — > ±v: a divergence which the e prescription does not 
cure. In contrast, the first term of (E3) does not give rise to a logarithm. For example, for 
m — 0, one gets 

(00\( lVz p + er 2 \00)= 1 - [ +1 d \ =1. (E5) 

IJ-i {tv z p + e) z p z 

We will now focus on the logarithmic divergence and so ignore the first term of (E3). The 
divergence is cut off only by 5C itself. That is, it is cut off when \iv z p\ becomes of order 
5c = 0(1), so that \v z \ ~ p _1 . The logarithmic divergence of (E3) is an artifact of our 
treatment of 5c as a perturbation. However, by proceeding with (E3) and treating v z as 
effectively cut off at 0(p _1 ), and so 5(v z p) as being smeared to have a width of A(i> 2 p) ~ 1, 
we may extract the coefficient of the logarithm. 

The log divergence can be seen by setting v z — v' z — and considering just the integral 
over azimuthal angles in (E3): 
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/ d<\) d<t>' ■ 
JO 



5c 2 (v,V) e im4 



f 2lT ,, . s CO: 

— 8 / d(A(f)) — 
Jo si 



cos 2 (A0) 
sin(A0)| 



im A<p 



(E6) 



There is then a log singularity associated with A0 — > 0. If one repeats the above for small 
(rather than zero) v z and v' z , one finds the dominant logarithmic behavior is 



/ d<f> d(f>' e- im<t> 5c 2 (v, v') e im * « -16 
Jo 



In 



v 2 - v' 



;-) m in 



1 



(E7) 



For even m, if we substitute back into (E3), treat the 5-functions as smeared over 
v z p ~ 1, and ignore non-logarithmic corrections, then the logarithms of (E7) becomes simply 
ln(l/p)'s, giving 



(o) ~ (2m + 1) ip-tO)? 2 In. - ( 2 ^ + 1 ) !! (2m-l)!! 2 
*£ m (p)~ (2m) , [P m (0)j ^lnp- ^ -^lnp. 



(E8) 



For odd m, we instead take the principal part terms of (E3), and we need to be more careful 
to separate the (cut-off) integral of the joint overall scale of v z and v' z from that of the 
relative scale (3 = v z /v' z : 



^ (2m + 1) 2 1 r+oo dv z r+- 

(2^)T [Pm (0)J V loo ~ loo dP ln 

(2m + 1)!! (2m- 1)!! 2 
* (2^! ^ lnp - 

Our final result is then 

(2m + 1)!! (2m- 1)!! 



V z + V' 



1 + p 



1-/9 



S TO (p) = 



(2m)! 



- +( _ r | lnp+0 i 



(E9) 



(E10) 



We do not know how to determine the non-logarithmic 0(p~ 2 ) piece analytically, but fits 
to numerical evaluation of the continued-fraction formula (2.116) at large p give the leading 
correction to (E10) to be lnp — > lnp + h m with h ~ —0.8, hi ~ 1.16, and h 2 — —4.3. 
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